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For  repairable  systems  with  standbys,  several  candidates  for 
measuring  reliability  are  (1)  the  time  to  system  failure,  (2)  the 
number  of  repairs  completed  before  failure  and  (3)  the  total  time 
spent  on  repair  before  failure.   In  this  research,  some  general 
procedures  are  developed  for  studying  these  and  other  reliability 
characteristics  of  systems  with  standbys.   Emphasis  is  placed  on  a 
class  of  repairable  systems  kno™  as  k-out-of-n  systems,  1  £  k  £  n, 
where  n  units  in  parallel  redundancy  are  serviced  by  a  single 
repairman  and  system  failure  occurs  when  k  units  are  simultaneously 
inoperable  for  the  first  time. 

Assuming  X.  and  F.(-)  are,  respectively ,  the  failure  rate  and 
general  repair  distribution  for  unit  i,  the  supplementary  variable 
technique  is  used  to  find  the  transform  of  the  time  to  system  failure 
distribution  for  the  k-out-of-n  system.   Using  the  principle  of 
regeneration,  an  alternative  derivation  is  also  obtained  for  the 
2-out-of-n  system. 


IX 


A  conditional  transform  approach  is  applied  to  the  2-out-of-n 
system  and  the  two-unit  standby  redundant  system.   For  each  system, 
transforms  of  distributions  are  obtained  for  the  time  to  system 
failure,  the  idle  time  of  the  repairman  and  the  time  spent  on  repair. 
In  the  case  of  exponential  failure  and  Erlangian  repair  capabilities, 
the  generating  function  for  the  distribution  of  the  number  of  renewals 
occurring  during  the  life  of  each  system  is  found. 

Numerical  results  are  presented  for  the  2-out-of-3  and 
3-out-of-4  systems  for  several  different  repair  disciplines.   Numerical 
methods  for  evaluating  moments  of  distributions  expressed  in  terms  of 
complex  integrals  are  also  discussed.   Although  specific  results  are 
not  derived,  the  value  of  stochastic  integrals  for  investigating 
attributes  of  redundant  systems  is  clearly  demonstrated. 


CHAPTER   1 
AN   OVERVIEVJ 

Introduction 

High  speed  computers,  sophisticated  weaponry  and  faster 
means  of  communication  and  transportation  have  focused  attention  in 
the  scientific  community  on  the  design  and  manufacture  of  systems 
which  are  both  proficient  and  reliable.   Most  engineers  would  agree, 
however,  that  unless  reliability  is  incorporated  in  the  design  phase, 
the  final  product  will  suffer  regardless  of  corrective  maintenance 
policy.   Thus,  it  is  the  ability  to  measure  reliability  that  concerns 
the  engineer  and  designer  of  today. 

For  repairable  systems  with  standbys,  several  candidates  for 
measuring  reliability  are  (1)  the  time  to  system  failure,  (2)  the 
number  of  repairs  completed  before  failure  and  (3)  the  total  time 
spent  on  repair  before  failure.   In  this  dissertation,  some  general 
procedures  are  developed  for  studying  these  and  other  reliability 
characteristics  of  systems  with  standbys.   Emphasis  is  placed  on  a 
class  cf  repairable  systems  known  as  k-out-of-n  systems,  1  <_  k  <_  n, 
whera  n  ur.its  (or  machines)  in  parallel  redundancy  are  serviced  by  a 
sin3l:;  repriirr.ian  and  system  failure  occurs  vjhen  k  units  are  simul- 
t-irieously  ii^operable  for  the  first  time.   The  often  referred  to  series 
ani  paralleJ  arrangements  are  the  special  cases  k  =  1  and  k  =  n, 
respectively.   Because  k-out-of-n  systems  (or  combinations  thereof) 
model  a  large  class  of  redundant  system.s,  Lhis  research  encompasses  a 
study  of  the  attributes  of  such  systems. 


Previous  Research  Results 

The  body  of  methods  and  ideas  used  in  studying  characteristics 
related  to  the  life  of  a  system  is  referred  to  as  mathematical 
reliability  theory.   Although  a  relatively  new  area  of  study,  reliability 
theory  has  already  been  applied  to  problems  involving  power  systems 
[1],  [2],  [3],  [4],  [5]  and  the  regulation  of  traffic  [6],  [7],  as 
well  as  to  multicomponent  structures  [8],  [9],  [10],  [11]  and  [12]. 

Takacs  [13]  in  1957  and  Morse  [14]  in  1958  discussed  properties 
of  a  system  where  several  machines  were  serviced  by  a  single  repairman, 
and  each  machine  was  subject  to  exponential  failure.   Limiting 
distributions  for  the  number  of  machines  working  were  derived  by  Takacs 
for  general  repair  and  by  Morse  for  exponential  repair. 

One  of  the  earliest  to  treat  first  passage  times  [15]  in  the 
context  of  reliability  theory  was  Gaver  [16]  in  I960.   For  a  birth- 
death  process  with  constant  transition  rates,  Gaver  derived  a  formula 
for  E[T.,],  j  >  i,  where  T..  is  the  first  passage  time  to  go  from 
state  i  to  state  j.   In  1962,  Belyayev  [17]  found  the  mean  time  to 
system  failure  (MTSF)  for  a  system  composed  of  two  units  in  parallel 
redundancy,  both  of  which  are  operating  initially.   Belyayev  assumed 
that  failures  were  exponential,  repair  x>7as  general,  and  the  system 
failed  when  both  units  were  inoperable  for  the  first  time.   In  1963, 
Gaver  [18]  solved  the  same  problem  considered  by  Belyayev  but  used   a 
modified  approach.   Both  Belyayev  and  Gaver  used  the  Cox  [19]  method 
of  supplementary  variables  to  analyze  the  process.   Later  in  1963, 
McGregor  [20]  derived  formulas  for  the  MTSF  for  the  more  general 
k-out-of-n  system  where  both  failures  and  repairs  were  exponential. 


Thiruvengadam  and  Jaiswal  [21]  in  1964  and  Jaiswal  [22]  in 
1968  used  discrete  transforms  to  investigate  the  k-out-of-n  system  with 
general  repair  capabilities.   In  addition  to  several  other  system 
attributes,  Thiruvengadam  and  Jaiswal  found  the  distribution  of  a  busy 
period  generated  by  i  units,  but  they  did  not  consider  properties  of 
the  time  to  system  failure. 

In  1966,  Muth  [23],  Htun  [24]  and  Srinivasan  [25]  discussed  a 
standby  redundant  system  composed  of  two  units.   Muth  found  the  OTSF 
when  the  repairman  has  the  capability  of  completing  n  repairs.   As 
n  -^  -,  Muth  showed  that  the  MTSF  agreed  with  equation  (20)  in  Gaver 
[18].   Htun  assumed  repair  and  failure  rates  were  constant  but  different 
for  each  unit  and  derived  the  MTSF.   Srinivasan  assumed  general  failure 
and  repair  capabilities  and  derived  the  transform  of  the  time  to 
failure  distribution. 

Also  in  1966,  Liebowitz  [26]  and  Downton  [27]  obtained  some  new 
results  for  parallel  redundant  systems.   Considering  the  two-unit 
redundant  system  treated  also  by  Belyayev  and  Gaver,  Liebowitz  plotted 
the  ratio  (called  the  improvement  factor) 

MTSF 
,-1 

/v 

for  various  repair  disciplines,  where  x"   is  the  expected  time  to 
failure  for  e.-ich  unit.   He  concluded  that  the  difference  between  MTSF 
with  repair  to  that  without  was  independent  of  the  repair  distribution 
itself  and  d-^ended  only  on  the  quantity  E[R] ,  where  E[R]  is  the  expected 
rcoair  uime  .^or  each  component.   Downton  used  the  properties  of  semi- 


Markov  processes  to  analyze  a  parallel  redundant  system  composed  of 

n  units.   Assuming  exponential  failure  and  general  repair,  Downton 

obtained  the  Laplace  transform  of  the  distribution  of  T,   ,  the  time 
^  k,n 

to  failure  for  a  k-out-of-n  system,  and  E[T    ]  for  n  >_  k,  k  =  2,3,4. 

ix. }  n        ' 

The  role  of  priorities  in  a  slightly  more  complex  system  than 
those  described  above  was  treated  by  Natarajan  [28]  in  1967.   Two 
paralleled  radars  working  in  conjunction  with  two  paralleled  computers 
describe   the  system  which  is  considered  operative  so  long  as  at  least 
one  computer  and  one  radar  are  working.   Natarajan  derived  the  distri- 
bution of  the  time  to  failure  and  the  MTSF  assuming  failures  and 
repairs  were  exponential.   Values  of  the  MTSF  were  compared  for  the 
case  of  no  priorities  and  when  the  preemptive-resume  discipline  was 
assumed.   In  his  doctoral  dissertation  [29]  completed  in  1968,  Natarajan 
treated  a  large  class  of  reliability  models.   Both  single  and  multiple 
repair  facilities  were  discussed  for  a  system  composed  of  a  basic  unit 
supported  by  n-1  standbys,  where  the  standbys  may  deteriorate  in 
storage  (sometimes   called  warm  standbys) . 

Also  in  1968,  Mine,  Osaki  and  Asakura  [30]  derived  the  MTSF  for 
a  2-out-of-n  system  by  considering  first  passage  times.   The 
improvement  factor  was  evaluated  and  its  asymptotic  behavior  was 
studied. 

In  1969,  Epstein  [31]  considered  a  two-unit  parallel  redundant 
system  where  failure  was  exponential  (but  with  a  different  failure  rate 
for  each  unit)  and  repair  was  general.   Transforms  of  distributions 
were  obtained  for  the  time  to  system  failure,  the  total  time  spent  on 
repair  and  the  free  time  of  the  repairman. 


In  1970,  Rao  and  Natarajan  [32]  presented  a  survey  paper  on 
the  reliability  of  systems  with  hot,  warm  and  cold  standbys  for  both 
parallel  redundant  systems  and  standby  redundant  systems.   Rao  and 
Natarajan  remarked  that  the  problem  of  finding  the  distribution  of  the 
time  to  system  failure  for  a  k-out-of-n  system  has  not  been  solved  when 
either  the  repair  or  the  failure  distribution  is  not  exponential. 
Also  in  1970,  Osaki,  Muth  and  Mazumdar  analyzed  some  two-unit  systems. 
Osaki  [33]  treated  a  two-unit  standby  redundant  system  and  found  the 
MTSF  using  a  state  transition  diagram  (or  signal  flow  graph)  approach. 
Osaki  [34]  applied  the  integral  equation  of  renewal  theory  to  four  two- 
unit  redundant  systems  (parallel  redundancy,  standby  redundancy  with 
and  without  priorities  and  standby  redundancy  with  noninstantaneous 
switchover).   The  transform  of  the  time  to  system  failure  was  derived 
in  each  case.   Muth  [35]  took  advantage  of  the  regenerative  properties 
(discussed  extensively  by  Smith  [36])  of  the  two-unit  parallel  redundant 
system,  and  derived  the  time  to  failure  distribution  and  the  MTSF  in 
an  elegant  manner.   Mazumdar  [37]  considered  a  two-unit  redundant 
system  where  the  operating  unit  and  the  standby  may  have  different 
failure  rates.   Mazumdar  found  the  MTSF  when  the  detection  of  failure 
for  the  standby  and  the  operating  unit  may  not  be  found  instantaneously. 

In  1971,  Gaver  and  Luckew  [38]  treated  what  they  termed  an 
accumulation  model,  where  a  system  experiences  accumulating  degrading, 
but  not  fatal,  failures.   Joint  transforms  were  derived  relating  total 
time  spent  on  repair  and  the  number  of  repairs  made.   Later  in  1971, 
Branson  and  Shah  [39]  considered  a  two-unit  parallel  redundant  system 
where  each  unit  had  a  different  but  constant  failure  rate,  and  repair 


distributions  were  different  but  general  for  each  unit.   Using 
properties  of  semi-Markov  processes,  Branson  and  Shah  dervied  the 
MTSF  for  the  two-unit  system  and  discussed  some  of  the  difficulties 
encountered  when  the  semi-Markov  model  was  applied  to  a  three-unit 
system  under  similar  assumptions. 

Problem  Statement 

Although  many  authors  discuss  properties  of  T,  the  time  to 
system  failure,  other  measures  of  system  reliability  have  been 
generally  neglected  in  the  literature.   This  dissertation  will  exhibit 
techniques  for  finding  transforms  of  distributions  for  the  time  spent 
on  repairs  during  T,  the  time  that  the  repairman  is  idle  during  T  and 
the  number  of  breakdowns  prior  to  system  failure.   In  addition,  assuming 
random  (exponential)  failure,  general  repair  capabilities  and  one 
repairman,  the  supplementary  variable  technique  of  Cox  [19] ,  Belyayev 
[17]  and  Gaver  [18]  will  be  extended  to  allow  analysis  of  the  k-out-of-n 
system.   A  comparison  will  be  made  between  the  supplementary  variable 
approach  used  in  this  research  and  the  semi-Markov  model  employed  by 
Doumton  [27]  and  Branson  and  Shah  [39]. 

Research  Objectives 

Specifically,  the  objectives  of  this  dissertation  can  be 
summarized  as  follows: 

(1)   For  the  k-out-of-n  system,  to  find  the  transform  of  the 
distribution  of  T  assuming  that  ).      and  F  (•)  are,  respectively,  the 
failure  rate  and  general  repair  distribution  for  machine  m,  1  <  m  <  n. 


(2)  For  the  n-out-of-n  system  with  X   =  X  and  F  (•)  =  F(0, 

■'mm 

1  <_  m  <_  n,  to  find  a  general  expression  for  the  MTSF. 

(3)  For  the  2-out-of-n  system  and  the  two-unit  standby 
redundant  system,  to  derive  transforms  of  distributions  for  the  time 
spent  on  repairs  during  T,  the  time  that  the  repairman  is  idle  during 
T  and  the  number  of  repairs  completed  prior  to  system  failure.   Also, 
assuming  exponential  failure  and  Erlangian  repair  capabilities,  to 
analyze  the  distribution  of  T  for  each  system  and  apply  the  results  to 
a  variation  of  the  accumulation  model. 

(4)  To  provide  the  reader  with  insight  into  solutions  of 
problems  which  are  as  yet  unsolved;  namely,  the  analysis  of  the  2-out- 
of-2  system  when  failures  and  repairs  follow  general  distributions  and 
the  application  of  stochastic  integrals  to  redundant  systems. 

Research  Profile 

In  Chapter  2,  the  transform  of  the  time  to  failure  distribution 
is  obtained  for  some  k-out-of-n  systems  and  a  comparison  is  made 
between  the  supplementary  variable  approach  and  the  semi-Markov  method. 

Chapter  3  concerns  the  2-out-of-n  system  and  the  two-unit 
standby  redundant  system.   A  conditional  transform  approach  is  applied 
and  the  distribution  of  the  time  to  system  failure  for  both  systems 
is  analyzed. 

Chapter  4  deals  with  applications  of  the  research. 

Conclusions  and  areas  for  future  research  are  discussed 
in  Chapter  5. 


CHAPTER  2 
THE  k-OUT~OF-n  SYSTEM 

Introduction 

In  this  chapter,  the  supplementary  variable  technique  is  used 

to  find  the  transform  of  the  time  to  system  failure  for  the  k-out-of-n 

system.   After  introducing  some  preliminaries,  the  k-out-of-n  system 

is  treated  first,  when  A   =  X  >  0,  F  (•)  =  F(0,  1  Ji  "i  £  n>  and  second, 

when  X      and  F  (•)  are  distinct  for  each  m.   A  comparison  between  the 
mm  ^ 

supplementary  variable  method  and  the  semi-Markov  approach  is  also 
discussed . 

Assumptions,  Definitions  and  Notation 

The  following  will  be  assumed: 

(1)  n  units  (machines)  are  in  parallel  redundancy  and  the 

time  to  failure  of  each  machine  is  independently  and  exponentially 

distributed.   The  failure  rate  of  machine  m  is  A   >  0,  1  <  m  <  n. 

m    '   —   ~ 

(2)  At  t  =  0,  all  units  are  operative. 

(3)  There  is  only  one  repairman  (with  unlimited  service 
capabilities)  and  failed  units  queue  up  for  service  on  a  "first-come, 
first-served"  basis.   Machine  m  is  serviced  according  to  a  general 
repair  distribution,  F  (•),    and  repair  times  are  independent  random 
varic.bi.es  vhich  are  also  independent  of  the  failure  times. 

i'i)      Repaired  units  are  put  back  into  operation  if  k-1  or 
less  units  are  in  a  failed  state.   When  k  units  are  simultaneously 
inoperable  for  the  first  time,  the  system  falls. 


The  following  notation  will  be  used: 

(a)  T,     -   a  continuous  random  variable  which  for  any 
k-out-of-n  system  represents  the  length  of  time  elapsed  before  k 
machines  are  simultaneously  inoperable  for  the  first  time,  assuming 
all  n  machines  are  running  at  t  =  0;  0  <  T,    <  <». 

(b)  R  .   -   continuous  random  variables  which  represent  the 
time  needed  to  repair  the  i-th  breakdown  of  machine  m;  0  <  R  .  <  °°, 
Tn=l,2,  .  .  .  ,n;  j=l,2 

(c)  y  -  a  continuous  random  variable  which  at  any  moment  t 
represents  the  time  already  spent  on  the  repair  of  the  m.achine  under- 
going service  at  that  moment;  0  <_  y  <_  t  <  °°. 

(d)  i   -   a  discrete  random  variable  which  at  any  moment  t 
represents  the  number  of  machines  broken  do'Am  at  that  moment; 

(e)  E  -   the  event  "T,    >  t." 

k,n 

(f)  Let  F  (•)  denote  the  common  distribution  function  for  the 

m 

random,  variables  R  .;  i.e.,  F  (r)  =P[R  ,  <  r],  m=l,2,...,n;  i=l,2,.... 

mj         m         mj  —        >  >    .  »  j   »  . 

(g)  Put  P    (t)  =  P[E]  and  hence,  P    (0)  =  1  by  assumption  (2). 

(h)   Let  W,   (s)  denote  the  Laplace  transform  of 
K.,n 


^         (t)  •  -i  e    W    (<i)    = 


e  ^y   (t)dt. 
k,n 


(i)   Let  h,   (t)  and  h,   (s)  represent  the  density  of  T, 

K.,n         K.,n  -^     k,n 

and  its  Laplace  transform,  respectively. 

(i)   Let  X(t)  represent  the  state  of  the  k-out-of-n  system  at 
anj'  nomen  t  t ;  e.g., 

(i)        X(t)    =    (i,z),    l<i<k<n,    z>0,   means    that   at 
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moment  t,  there  are  i  machines  awaiting  repair  (including  the  machine 
presently  being  serviced)  and  an  amount  of  time  z,  0  <_  z  <_  t,  has 
already  been  spent  on  the  repair  of  the  machine  presently  being  serviced, 

(ii)   X(t)  =  (0,0)  means  that  at  moment  t,  there  are  no 
machines  broken  down  (i.e.,  no  machines  awaiting  or  being  repaired)  and, 
consequently,  no  time  has  been  spent  on  the  repair  of  a  machine. 

(k)   Put  pQ(t)  =  P[X(t)  =  (0,0)];  hence,  p^(t   =  0)  =  1  by 
assumption  (2) . 

(1)   Let  p_j^(t,z)  represent  continuous  and  dif f erentiable 
functions  such  that  p.(t,z)dz  =  P[X(t)  =  (i,y)],  1  1  i  £  k  <^  n, 
z  <  y  <^  z  +  dz. 

(m)   Define 


and 


Vo(s)  = 


V.(s) 


~st   ,  ,  , 
e   PQ(t)dt 


-St, 

e    ( 


p^(t,z)dz)dt,   1  £  i  ^  k-1 


(n)   By  definitions  (k) ,  (1)  and  (g) 


^,n(^)  =  PO^^)  +  T 


p^(t,z)dz 


and  hence,  using  (m)  and  (h) 


k-1 


\   n^^^  =  I        ^•(^)'   1  <  k  <  n 


11 


The  Case  ofX   =AandF(0=F(-),llniin 

za in 

Preliminaries 

From  the  above  assumptions  and  definitions,  transition 
probabilities  for  X(-)  during  the  time  interval  (t,  t  +  dt)  depend 
only  on  the  state  of  the  system  at  moment  t.   In  particular,  when 

A   =  A  >  0     and     F  (•)  =  F(-),  1  1  m  1  n         (2.1) 
m  ni 

the  transition  probabilities  can  be  described  as  follows: 
(i)    P[X(t  +  dt)  =  (l,0)|x(t)  =  (0,0)]  =  nXdt  +  o(dt) 
(ii)    P[X(t  +  dt)  =  (0,0)!x(t)  =  (0,0)]  =  1  -  n>dt  +  o(dt) 
(iii)   P[X(t  +  dt)  =  (j  +  1,  z  +  dt)lx(t)  =  (j,z)] 

=  (n  -  j)Xdt{(l  -  F(z  +  dt))/(l  -  F(z))}  +  o(dt),  1  1  j  lk-1 
(iv)    P[X(t  +  dt)  =  (j,z  +  dt)|x(t)  =  (j.z)] 

=  (1  -  (n-j)Xdt){(l  -  F(z  +  dt))/(l  -  F(z))}  +  o(dt),  1  1  j  1  '-c-l 
(v)    PLX(t  +  dt)  =  (j-l,0)|x(t)  =  (j,z)] 

=  (1  -  (n-j)Xdt){(F(z  +  dt)  -  F(z))/(1  -  F(z))}  +  o(dt),  IfJlk-l 

The  2-"0ut-cf-n  System 

Consider  the  2-out-cf-n  system.   For  n=2,  a  typical  realization 
of  the  process  XCt)  in  terms  of  the  random  variables  y  and  i  is 
presented  in  Figure  2.1. 
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2   -■ 

1   -- 
0 


'l       '^l-'^l  ^2    ^^3 

Figure  2.1.   Sample  function  of  the  process  X(t)  =  (i,y), 
i  =  0,  1,  2,  y  >_  0.   t.  is  the  moment  at  which  the  j-th 

breakdo^^m   occurs  and  r^  is  the  time  needed  to  complete  the 

first  repair.   The  second  repair  is  not  completed  before 
the  absorbing  state  (2,z)  is  reached,  and  thus,  T    =  t„. 

Two  quantities  of  interest  are  the  distribution  of  the  random 

variable  T^   and  its  expectation,  E[T_   ].   In  this  regard,  from 
2,n  "  2,n  "    ' 

definitions  (g)  and  (h)  it  is  seen  that 


^2.n<^)  =  dl  (1 


P,   (t))  = 
2,n 


dt   2,n 


(2.2) 


and  since  P    (t=0)  =  1,  from  (2.2) 
z  ,n 


h    (s)  =  1  -  sW„  ^(s) 
2,n  z,n 


(2.3) 


Furthermore,  assuming  E[T„   ]  <  °°,  from  (2.3) 

z ,  n 


E[T  „]  =  (1  -  h_  ^(s))/s 
2,n  2,n 


s=0 


=  W    (s) 
Z,  n 


s=0 


(2.4) 


To  find  W   (s) ,  first  consider  the  state  transition  equations. 
/ ,  n 

Using  definitions  (k)  and  (1),  transition  probabilities  (i)-(v)  and  the 
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fact  that  (2,z),  z  ^0,  are  absorbing  states,  the  follov;ing  relation- 
ships hold  for  n  _>_  2: 


p^Ct  +  dt)  =  PQ(t)(l  -  nXdt)  + 


+  [  p^(t,z)(l  -  (n  -  DXdt)  ^^^  1  -^F(z)^^^^  ^^  "^  °^^'^ 


0 


(2.5a) 


p  (t  +  dt,z  +  dt)  =  p^(t,z)(l  -  (n 


p  (t  +  dt,0)dt  =  p  (t)nXdt  +  o(dt) 


(2.5b) 
(2.5c) 


Po(0)  =  1 


(2.5d) 


Defining  g  (t,z)  by  the  relationship 


p^(t,z)  =  g^(t,z)e  ^"^   ^^^'^  (1  -  F(z)) 


(2.6) 


and  using  the  fact   that  for  g(")  continuous  and  F(*)  absolutely 
continuous 


lim 
h->0 


r    ^    F(x  +  h)  -  F(x)  , 
g(x) dx 


J 


g(x)dF(x) 


(2.7) 


0  0 

equations  (2.5a)  -  (2.5d)  may  be  rewritten  as 


This  can  be  shown  by  using  the  properties  of  continuous  and 
absolutely  continuous  functions  and  appealing  to  the  Lebesgue  Dominated 
Convergence  Theorem. 


^^0 


e  ^"  ^^^^  g^(t,z)dF(z),pQ(0)  =  1 


(2.8a) 


'^1   ^^1 

— —  +  T— ^  =  0 

3t     dz 


p^Ct.O)  -  nXp^Ct) 


(2.8b) 
(2.8c) 


To  obtain  (2.8a),  substitute  (2.6)  into  (2.5a),  use  (2.7) 
and  take  the  limit  as  dt  ->  0.   To  obtain  (2,8b),  first  substitute  (2.6) 

into  (2.5b).   Noting  that  q-(^-^^^<^^   =  1  -  (n  -  l),\dt  +  o(dt)  obtain 
g  (t  +  dt,  z  +  dt)  =  g  (t,z)  +  o(dt).   Adding  in  and  subtracting  out 
the  quantities  g  (t,z  +  dt),  g-,  (t  +  dt,z)  and  g  (t,z)  in  the  preceding 
equation,  dividing  by  dt  and  taking  the  limit  as  dt  -*■  0  yields  (2.8b). 

Equations  (2.8a),  (2.8b)  and  (2.8c)  will  now  be  used  to  find 
W    (s) .   Before  doing  so,  note  that  by  definition  (n) 

^  }  XI 


W„      (s)    =  V^(s)   +  V    rs),      n  >   2 
/ ,  n  U  1  — 


(2.9) 


e   ^^dF(z) 


Also,    let   the   l.aplace-Stieltjes   transform  of   F(z)    be   <^is) 
and   define 


ij;   (s)    =   (J(s  +  jA)-,      j    >_  1 


Nov;,  W^   (s)  can  be  found  from  equations  (2.8a),  (2.8b)  and 
.1,  n 

(2.8c.)  as  follows.   From  (2.8b) 


g^(t,z)  =  gj^(t  -  z) 


(2.10) 


Combining  (2.6)  and  (2.10)  with  (2.8c) 
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g^(t)  -  riXp^Ct) 


(2.11) 


Defining  G  (s) 


—St 

e   g  (t)dt,  take  the  Laplace  transform  of  (2.8a) 


0 


and  (2.11)  and  obtain  respectively 


1  +  iij    (s)  •  G  (s) 

V  (3)  =  ^i 1 

0^"^        s  +  nX 


and 


or 


G^(s)  =  nXVQ(s) 


V  (s)  =  {s  +  nX(l  -  i^      As))} 
U  n-i 


-1 


(2.12) 


From  (2.6)  and  (2.10) 


V^(s)  =  I   e  ^^  I  e"^''  ^^^^g^(t  -  2)(1  -  F(z))dz)dt 


0 


0 


G,(s) 
s  +  (n-l)X  (^  -  Vl^^^) 


nX(l  -  ^      As)) 
n— i 


(s  +  (n-l)X){s  +  nX(l  -  i>      ,(s))} 

n-X 


(2.13) 


Thus,  from  (2.2),  (2.3),  (2.4),  (2.12)  and  (2.13)  it  follows 
tliit  for  the  2-cat-of-n  system 

_-,  nX(l  -  4)  _  (s)) 

W„  ^(s)  -  {s  +  n\(l  -  -.J;   .(s))}  ^  +  .      ^    ,      ..... — ",  ,,     t-tyT 

-^.n  n-1  (s  +  (n-j.)X){s  +  nX(l  -  ^      .i.s))) 

n-i 


(2.14) 


h..  ^(s)  =  1  -  sW.   (s) 
z, n  z, n 


and 
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E[T„   ]  =  W-   (0)  =  7-^T--  +   WT  '*'  ■  ■  ■  -s 
2,n     2,n      (n-l)A   nA(l  -  <p      ,) 

n-i 


(2.15) 


where  <!:_.=  $((n-l)>.).   Relation  (2.15)  agrees  with  equation  (5.19a) 
in  Down ten  [27] . 

The  3-Out-of-n  System 

For  the  3-out-of-n  system,  U   (s)  can  be  found  in  a  similar 

3,n 

manner.   From  definition  (n) 


W-  ^(s)  =  I     V.(s) 
^'"      i=0  ^ 


(2.16) 


and  as  before 


iad 


h„  ^(s)  =  1  -  sW-   (s) 
->,n  J,n 


E[T    ]  =  W    (0) 
J  ,n     J  ,n 


Defining  g,(t,z)  and  q„(t,z)  by  the  relationships 


p^(t,z)  =  e  ^"^  ^^^^  g;L(t,z)(l  -  F(z)) 


(2.17) 


q2(t.z)  =  P2(t,z)/(1  -  F(z)) 


(2.18) 


the  state  transition  equations  for  the  3-out-of-n  system  may  be 
rewritten  as 


dT  =  -^^0^^^  -^ 


e  ^"^  ^^''^  g^(t,z)dF(z),pQ(0)  =  1 


(2.19a) 


9g-,    3gi 
~  +  _~-  =  0 

d  t      dZ 


(2.19b) 
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T-^  +  ^-=  -(n-2)Aq,(t,z)  +  (n-DXe  ^"""^^^^gAt  ,z)  (2.19c) 


p^(t,0)  =  nXp^Ct)  + 


q2(t,z)dF(z)  (2.19d) 


P2(t,0)  =  0  (2.19e) 

Again,  equation  (2.7)  was  used  in  obtaining  equations  (2.19a)  and 

(2.i9d).   From  (2.19a)  -  (2.19e),  W    (s)  can  be  obtained  as  follows. 

J ,  n 


From  (2.19b) 


g^(t,z)  =  g^(t  -  z)  (2.20) 


Substituting  (2.20)  into  (2.19c)  and  solving  yields 

q2(t,z)  -.   e-(^-2)^2  g^(t  -  z)  -  (n-l)e"^"-^^^^g^(t  -  z)      (2.21) 

where  g„(*)  is  an  arbitrary  function  of  t  ~  z.   Equations  (2.18)  and 
(2.19e)  imply  that  q_(t,0)  =  0,  and  thus,  (2.21)  becomes 

/_   \    /   i\   /_    \f    -(n-2)Xz    -(n-l).\z-,  ,^    „,s 

q2(t,z)  =  (n-l)g^(t  -  z){e         -  e  ^    '   }  (2.22) 

Using  (2.17)  and  (2.20) 


p^(t,0)  =  g^(t,0)  =  g^(t)  (2.23) 

Combining  (2.22)  and  (2.23)  with  (2.19d)  yields 

t 

/^\  ^   /■^\     1  /   IN     /^     \r    -(n-2)Xz     -  (n-1)  Az-i  ,  ,w  n 

g^it)    =  nApQ(t)  +  (a-l>   g^(t  -  z){e         -  e  ^      }di'(z) 

0 

(2.24) 

V  (s)  m.ust  be  found  from  equations  (2.19a)  and  (2.24).   To  this  end, 

after  taking  the  Laplace  transform  of  (2.19a)  and  (2.24;,  obtain  respectively 
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1  +  G  (s)  •  ^        (s) 
V  (s)  = i -^^ 


(2.25) 


and 


hXVqCs) 
G^(s)  =  YT:— n-l)i^^    rs)  +  (n-l)^„  ,  G) 


'n-2 


•^n-l 


(2.26) 


where  again 


G^(s)  = 


-St    ,    ■,  , 
e     gj^U)dt 


and 


1^.  (s)  = 


e  ^^^^^^^dF(t)  =  *(s  +  jX),   j  >  1 


From  (2.25)  and  (2.26) 


1  -  (n-2)^   „(s)  +  (n-l)(!;   .  (s)  ■ 
,.  /  N  ^  r}-2. n-1 

O^^'^    (s  +  nX){l  -  (n-l)jj  .(s)  +  (n-D^j   ^(s)}  -•  n\^      .  (s) 


n-z 


'n-1 


'n-1 


(2.27) 


A].so,  from  (2.17)  and  the  definition  of  V^  (s) 


V^(s)  =1   e  °^ 


e  ^"  ^^^^  g^(t  -  z)(l  -  F(z))dz)dt 


G/s) 
s  -^  (n-l)X  (1  -  Vl(  =  )> 


(2.28) 


Similarly,  from  (2.18),  (2.22)  and  the  definition  of  V  (s) 


r 

V_(3)  =  (n-1)  !  e 
2  J 

0 


(n-i;G^(s) 


-St 


/I    /     \r  -(n-2)Xz    -(n-l)Xz, /.    „,  .,,  ..^ 
(  I  g^(t  -  z){e         -  e  ^   ^   }(1  -  F(z))dz)dt 


0 


^  1  -  V2^  =  )    1  -  ^'^n-l^^^ 


(2.29) 


s  +  (n-2)X    3  +  (n-l)X 
Using  (2.27),  (2.?, 8),  (2.29)  and  (2.16).  one  finds  after  simplification 
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W^   (s)  =  [(s  +  nX)(n-l)(n-2)A  • 
j,n 

•  {X  -  (s  +  (n-l)A)C>   .(s)  +  (s  +  (n-2)>.)^^.   ,(s)}]/ 

n-/  n-i 

/[(s  +  nA){l  -  (n-l)i^   „(s)  +  (n-1)';;  As)}    -  nX^lJ      As)]  (2.30) 

n-/  n— 1  n— X 

and 

^3,n(^>  =  1  -  ^"3.n(^) 

1.1.      1  -  *n-l  +  V2 


E[To  „]  =  W^  _(0)  =  ,       ,.,    +   TT-^TTT  + 


■3,n^     3,n^''^    (n-l)X    (n-2)X  "  nA{l  -  (n-l)(^   .  +  (n-2)$  A 

n— /        n— i 

(2.31) 

where  ({> .  =  ij;.(0)  =  4>(j^)>  j  .^1-   Equation  (2.31)  agrees  with 
equation  (5.19b)  in  Downton  [27]. 

It  should  nov7  be  clear,  however,  that  although  E[T,   ]  can  be 

K.  J  il 

found  for  any  fixed  k,  a  general  expression  for  E[T    ]  cannot   be 
obtained  by  this  approach;   the  reason  of  course  being   that  the 
set  of  differential-integral  equations  which  result  from  the  state 
transition  equations  is  different  for  each  value  of  k.   A  similar 
difficulty  encountered  by  Downton  will  be  discussed  below. 

The  n-Out-of-n  System 

Although  unsuccessful  in  obtaining  a  general  expression  for 

E[T    ]  it  will  now  be  shown  that  the  problem  of  finding  E[T    ]  can 
K. ,  n  n ,  n 

be  reduced  to  solving  a  system  of  (n-1)  x  (n-1)  linear  equations. 
First,  for  n  >_  4  (the  n=2  and  n=3  cases  follow  from  equations  (2.15) 
and  (2.31)  above),  the  state  equations  are 
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P|^(t  +  (it)    =   PQ(t)(l  -  nXdt)   +    I      p^(t,z)(l   -    (n-l)Adt) 


0 
F(z  +  dt)    -   F(z) 


1   -   F(z) 


dz  +  o(dt)  (2.32a) 


p^(t  +  dt,    z  +  dt)    =  p^(t,z)(l   -    (n-l)Adt)    ^   1  !^F(t)'~  ^  °'^^^^ 


(2.32b) 


Pj^(t  +  dt,z  +  dt)    =   Pj^(t,z)(l   -    (n-k).\dt)    ^   I  ^-^FCz)^~  "^ 

+  p,     ^  (t,z)(n  -  k  +  DAdt   -^   :  ^^w^,'^"^    +  o(dt),    2    <   k   <   n-1 
k-1  1   -   FCz)  —       — 


(2.32c) 


p^(t  +  dt,z  +  dt)    =   p^(t,z)    +  p^_^(t,z)Xdt    ^    I   !^F(z)'^''''    ^   ''^'^^^ 

(2.32d) 

p^(t  +  dt,0)dt   =   pQ(t)nXdt  +    I      P2(t,z)(l  -    (n-2)Adt)    • 

0 

•   HL-±4^I_Z^)_  dz  -.  o(dt)  (2.32e) 

t 
p,^(t  +  dt,0)dt  =    f      pj^_^^(t,z)(l   -    (n-k-l)Xdt)    ^^^  t  -%^r^  "^^  " 
0 
+  o(dt),        2   <  k  <   n-2  (2.32f) 


p.(t  +  dt,0)dt   =   o(dt),      j=n-l,n  (2.32g) 


t-i 


,(0)    =   1  .  (2.32h) 


Proceeding   sin-.ilarly   as   above,    define   g    (t,z)    and    {q  (t,z)}             by 
tha   relationships 

p^(t:,z)    =   g^(t,z)e"^"~^^''^l   -  F(z))  (2.33) 

p    (t,z)    =  a    (t,z)(l   -  F(z)),      2  ^  k  <_  n-1  (2.34) 


Now,  using  again  the  fact  that  for  g(')  continuous  and  F(')  absolutely 
continuous 
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lim 
h^O 


,    -.    F(x  +  h)  -  F(x)  ,     [    (    \A'ci    \ 
g(x)  — ^^ ~ dx  =    g(x)dF(x) 


0  0 

equations  (2,32a)  -  (2.32h)  may  be  rev/ritten  with  the  aid  of  (2.33) 
and  (2.34)  as 

t 


e  ^^   ^^^^  g^(t,z)dF(z),pQ(0)  =  1 


0 


T—  +  T-"*^  =  0 

3 1    3z 


(2.35a) 


(2,35b) 


~-  +  T-^=  -  (n-2)Xq,(t,z)  +  (n-l)Xe  ^'^   ^^^^  g,(t,z) 

(J  t      dZ  Z  X 


(2.35c) 


--—  +  T—  =  -  (n-k)Aq,  (t,z)  +  (n  -  k  +  l)Xq,  .(t,z),  3  <  k  <  n-1 
ot     oZ  k.  K.-1  —   — 

(2.35d) 


3?    3p 

•^  +  7r  =  ^'In  i(t,z)(l  -  F(z)) 
o  t     a  Z       n- i 


p^(t,0)  =  nXpQ(t)  + 


q2(t,z)dF(z) 


(2.35e) 


?j^(t,0)  =     qj^^^(t,z)dF(z),   2  £  k  <_  n-2 

0 
p,(t,0)  =  0,   j  =  n-l,n 


(2.35f) 
(2.35g) 


Elquations  (2,35a)  -  (2,35g)  will  now  be  solved  and  the  results 

;ed  to  find  W   (s) ,   First,  let  the  solution  of 
n,  n 


^^^k  ,  ^\ 


—  +  — -  +  aq,  (t,z)  =  f(t,z),   2  ^  k  ^  n-1 
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be  of  the  form 


qj^Ct.z)  =  e  °^gj^(t  -  z)  +  q^(t,z) 


where  the  g,  (t  -  z)  are  functions  whose  existence  and  uniqueness 
follow  from  the  existence  and  uniqueness  of  the  {p,  (t,z)}   „ 

K.  K.  £m 

and  equation  (2.34),  and  q  (t,z)  is  the  appropriate  particular  solution, 
Theorem  2.1.   The  solution  of  (2.35c)  and  (2.35d)  in  terms  of  the 

Ui^ct  -  z)};;:^  is 


qj^(t,z)    =      I       (-1)^ 


k-1 

I 

i=0 


n-k+i 


-(n-k+i)Az  ,  ,  ,  , 

%-i^      ~   ^•^  (2.36) 


for   2  £  k   <   n-1,    n   >    4. 


Proof.   The  proof  consists  of  substituting  (2.36)  into  (2.35c)  and 

(2.35d)  directly. || 

Now,  from  (2.35e)  -  (2.35g),  (2.33)  and  (2.34),  it  follows 

that 

t 
g^Ct)  =  P3^(t,0)  =  nApQ(t)  +    q2(t,z)dF(z)  (2.37) 


q^._^(t,0)  =  pj,_^(t,0) 


0 


qj^(t,z)dF(z),   3  £  k  £  n-1 


(2.38) 


Vl^^'°^  =  Pn-l<^t'°)  =  0 


(2.39) 


Substituting  (2.36)  into  (2.37),  (2.38)  and  (2.39),  and  taking 
Laplace  transforms  of  the  resultant  equations  and  (2.35a)  yields  a 
system  of  n-1  equations  in  terms  of  the  n-1  unknowns,  V  (s) ,  G  (s) , 
02(5),. ..,G  _2(s),  n  >_  4,  where 


e  ^^g.(t)dt,  1  1  j  1  n-2 


0  °° 


G.(s)    = 


Remembering  that  'J^.(s)  =  |   exp{-  (s  +  jA)t}dF(t),  the  set  of   n-1 
equations  in  n-1  unknoi<ms  can  be  written  in  matrix  notation  as 


D  (s)  •  g  (s)  =  e, 

n       n      — X 


(2.40) 


where 


g,(s)' 


(Vq(s).  G^(s).  G^Cs), 


'V2^^)) 


e^  =  (1,0,0, ...,0) 


and,  dropping  the  s's  from  the  ii.'s,    reference  the  matrix  on  page  24. 

Thus,  for  every  fixed  n  >_   4,  equation  (2.40)  may  be  solved 
for  V  (s)  ,  G  (s)  ,.  .  .,G  _2(s)  in  terms  of  the  ^.  (s),  1  <_  j  1  n-1. 
Using  (2.36)  and  the  definitions  of  {V  (s)} 


G^(s) 

^1^^-^  "  s  +  (n-1) A 


(1  -  *„  ,(s)) 
n-1 


k-1 
V,  (s)  =  I      (-1)^ 
^      i=0 


n-k+i 
i 


Vi^^^^l-  Vk4-i(^)> 


n-1 
i=2 


'n-1 


i(-l)^G   .(s) 
n-1 


s  +  (n-k+i) X 

f  1  -  '^..(s)    1  -  ^,  (s) 


,  2  <  k  <  n-2 


s  +  X 


s  +  iX 


Now,  by  the  definition  of  W   (s)  and  the  above 

n,n 

k=2 


G  (s)(l  -  i'      ^(s)) 
^0^"^  ^     3  +  (n-l)X     ^ 


CS| 


24 


CO 

-a- 

^ .^ 

+ 

oo  o 


I 


PI 
I 

c 


I 


I 


CM  »!■ 
1   1 

a  a 

M 

a 

■9- 

CH 

•         CN 

sr 

^= 

c 

e 

1 

r-t 

. 

-3- 

— ^ 

i—^ 

CM 

1 

i-M  cn 

a 

c 

1  t 

G    G 

C4 

■  ^ . 

1 

<N 

c 

1 

.'-N 

a 

I 

a 


I 

G 

i 


I    o 


I 


a 


I 


iH 

a 

-3- 

r-l 

rH  n 

c 

1  1 

-9- 

c  c 

1 

iH 

+ 

-*• 

'  ' 

.-1  VT 

y' -^ 

1   1 

r-i    CN 

c  ts 

1    1 

a  c 

fH 
1 

r-i 

c 

/-v 

G 

C9 
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n-2 

k-1 

'n-k+i' 

+     I 

I 

i 

k=2 

i=0 

■ 

G^_.(s)(l-   Vk+i^^^^ 


n-1 

r 
+        I 

i=2 


i(-l)^G      .(s) 
n-i 


+    (n-k+i)X 


1  -   ^^(s)        1   -   '1'.  (s) 


+ 


s  +  A 


s   +   iX 


I 


)r,    after   simplification 


n-2 
W     _(s)    =  V^(s)    +      I 

j=l 


n,n 


0 


(-1)''"^    G.(s)    • 
J 


(n-j)(l  -  <:.^(s))      1  -  Vj(s> 

s  +  X  s  +    (n-j)X 


(2.41) 


where,  for  each  value  of  n  1  4.  Vq(s),  G^(s) , . . . .G^.^Cs)  are  computed 
from  (2.40).   Hence,  from  (2.40)  and  (2.41)  one  may  compute 


h   (s)  =  1  -  sW   (s) 
n,n  n,n 


and 


E[T    ]  =  W    (0) 
"■  n,n     n,n 


for  the  n-out-of-n  system. 

Taking  n  =  4,  for  example,  (2.40)  becomes 


fs  +  4X 


-'^3(3) 


-(1+3^1.3(5)) 


1^2(3) 


I 


0         3a-^^i^)+^^(s))     -a-^-i{s)^2M^(s)) 


fVQ(s) 


G^(s) 


G2(s) 


0 


(2.42) 


Solving  (2.42)  for  Vq(s),  G^(s),  6^(3)  and  evaluating 
vie  Ids 


them  at  s  =  0 
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(1  +  $  )(1  -  2-,  +  U   )    -   3?_(1  -  *,  +  ^_) 

V  CO)  = = I         i ,„  ,^. 

0^  ^    4X{(.l  +  2*3)  (1  -  2$^  +  2^2)  -  3(^2(1  -  <f^  +  c^^)}     ^2'^^^ 

1  -  2$^  +  24,2 
S^°^  =  Tl  +  2*3)  (1  -  2<^^  +  2^2>  -  3*2(1  -  *i  +  *3>         ^^■^^'' 

3(1  -  ^i  +  ^3) 
^2^°^  ^  a  +  2$-)(l  -  2*  +  2$,)  -  3$  ri  -  ^^  +  *  )         (2.45) 


where  *^  =  -li^^CO)  =  ^(iA),  i  ^  1.   Evaluating  W    (s)  from  (2.41) 
at  s  =  0  and  substituting  (2.43)  -  (2.45),  by  algebra 


^[^.4^  =  ^.4^°) 


(2.46) 


25  -  32^^  -  49^2  +  ^9*3  +  63*^^2  +  15<^2*3  ~  ^^*  *3 
~  12a{1  -  2*^  -  <^2  "*■  2*3  -  4*i<^3  +  *2*3  "^  3^^<?2^ 

Equation  (2.46)  agrees  with  equation  (5.19c)  in  Downton  [27]. 


The  Case  of  A  ,  F  (•)  Distinct,  1  <  m  <  n 
m   m —   — . 

Preliminaries 

When  the  restriction  imposed  by  equation  (2.1)  is  relaxed  and 

A   and  F  (•)  are  different  for  each  machine,  let  T,    be  the  time  to 
la      m  K  ,n 

system  failure.   It  will  be  shown  that  the  problem  of  finding 

Ere"^^k,n]  .=  ^   (g) 

and  the  MTSF,  denoted  now  by  E[T   ],  remains  solvable  for  fixed  values 
ot    k,  k  <_  n.   The  supplementary  variable  technique  is  still  applicable 
but  a  redefinition  of  states  is  necessary. 
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In  order  to  avoid  new,  or  at  best  cumbersome,  notation,  the 
subscript  i  and  the  functions  p^Ct),  p.(t,z),  g.(t,z),  V.(s),  G.(s), 
^.(s)    and  *.  will  continue  to  be  used,  even  though  they  may  carry 
slightly  different  meanings  than  in  previous  sections.   The  reader  is 
:tioned  to  make  note  of  this  whenever  these  quantities  are  defined 


cau 


below. 

The  2-Out-of-n  System 

For  this  more  general  2-out-of-n  system,  let  X(t)  be  the  state 
of  the  system  at  any  moment  t;  e.g., 

(i)    X(t)  =  (i,z),  1  <_  i  1  n.  z  >_  0,  means  that  at  moment  t, 
achine  num.ber  i  is  down,  and  an  amount  of  time  z,  0  £  z  ^  t,  has 
already  been  spent  on  the  repair  of  machine  i. 

(ii)   X(t)  =  (0,0)  means,  as  before,  no  machines  are  down. 

Letting 

p^Ct)  =  P[X(t)  =  (0,0)] 


m 


Pi 


(t,z)dz  =  P[X(t)  -    (i,y)],  1  1  i  1  n,  z  <  y  1  z  +  dz 


a=  I      X..   X^"^  =  I     X    1.1  i  in 
the  state  transition  equations  can  be  written  as  follows 


^(t  +  dt)  =  Po(t)(l  -  adt)  + 


.t  ^  ,    F.(z  +  dt)  -  F  (z) 

+   I    I    P,(t,2)(l  -  ^-  'it)  —l"^'.Tzl 


""  '  ,(i),,^^)  111" i::___l-'-  dz  +  o(dt)   (2.47a) 


p^(t  +  dt,Z  +  dt)  =  p^(t,2)( 


^,     1  -  F,.  (z  +  dt) 


1  -  a'  'dt) 


1  -  F.(z) 


+  o(dt)  ,  1  <  i  <  n 


(2.47b) 


p^(t  +  dt,0)dt  =  pQ(t)A^dt  +  o(dt),   1  1  i  1  n 


(2.47c) 


PqCO)  =  1 


(2.47d) 


)efining  {g.(t,2)}._   by  the  relationships 


_j(i) 
p^(t,z)  =  e  '    ^  g.(t,z)(l  -  F^(z)),  1  _<  i  1  n 


(2.48) 


equations  (2.47a)  -  (2.47d)  can  be  rewritten  with  the  aid  of  relation 
(2.7)  as 


__  =  -apQ(t)  +   I   J   e       g.(t,z)dF.(z),pQ(0)  =  1 
^"■'"  0 


(2.49a) 


3gi    3g. 

T—  +  r-i  =  0,   1  <  i  <  n 

o  t      dZ  —    — 


;j^(t.O)  =  ?.(t,0)  =  X^Py(t),   1  1  i  1  n 


(2.49b) 
(2.49c) 


For  1  <^  i  <_  n,  put 


Vo(s) 


e  ^''pQ(t)dt 


V.(s) 


-St, 

e    ( 


p.  (t,z)dz)dt 


G^(s)  =  j   e    g^(t)dt 


!^^(s) 


e  ^^dF.(t) 
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^^(s)  =--  J^(3  +  A^^^) 
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and 


n 


W 


2  n^^^  =   ^  ^1^^^ 
^'^      i=0 


Equations  (2.A9a)  -  (2.49c)  can  be  used  to  find  W   (s) ,  and  hence 

^ ,  n 

E[T„   ]  =  W„   (0),  as  follows.   From  (2.49b) 
z,  n     2 ,n 


g^(t,z)  =  g^(t  -  z),  1  1  i  1  n 


(2.50) 


Combining  (2.50)  with  (2.49c)  and  taking  Laplace  transforms  yields 


G^(s)  =  A^.Vq(s),  1  1  i  1  n 


(2.51) 


From  (2.50),  (2.51)  and  (2.49a),  it  follows  that 


Vq(s)  =  {s  +  a  -  I     \^i:^{s)} 
i-1 


-1 


(2.52) 


Using  the  definitions  of  V.(s),  (2.48),  (2.51)  and  (2.52) 


o=[ 


V.(s)  =  !   e  ^^( 


p.  (t,z)dz)dt 


»lV  (s) 

O  T  A 


X.d  -  ij;.(s)) 


(s  +  x''^^){s  +  a  ~   y   XJ;.(s)} 

i=l   ^  ^ 


,  1  <  i  <  n 


(2.53) 


By  (2.5  2)  and  (2,53) 


W„   (s)  =   y   V.(s) 


^.n 


i-0 


n 


=  {s  +  ct  -  y  x.tj;. (s)} 
1=1  ^ " 


_T  I     n  X  (1  -  ^..  (s)) 

^  1  +  y  -i       ^ 


i-l   s  +  X 


..  .(i) 


(2.54) 


~* 

and  hence   the  transform  of  the  density  of  T„   ,  h„   (s) ,  and 

2,n   2,n 

E[T„   ]  can  be  computed  from  (2.54)  as 
2,n 


h„  ^(s)  =  1  -  sW_  ^(s) 
2, n  /  ,n 


and 


(2.55) 
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E[T   J  =  W    (0) 
z,n  2  ,n 


=  {a  -   y  X.<i>.} 
i=l  ^  " 


-1 


n   X.  (1  -  <|).) 
i=l     X^^^ 


(2.56) 


where  (}> .  =  <^ .  (X        )  ,    i  ^  1.   As  a  check,  when  X.  =  X  and  F.(')  =  F(-), 
<ti^(s)  =  <t>(s),    X^^^    =  (n-l)X,  V^.(s)  =  (t>(s  +  (n-l)X),  a  =  nX  and  (2.54) 
becomes 


W    (s)  =  {s  +  nX[l  -  ^(s  +  (n-l)X)]} 
z,n 


-1 


1  + 


n  [1  -  (^(s  +  (n-l)X)] 
s  +  (n-l)X 


+ 


nX[l  -  <^(s  +    (n-l)X)] 


s  +  nX[l  -  (i)(s  +  (n-l)X)]    (s  +  (n-l)X){s  +  nX[l  -  <J  (5+(n-l)>,)  ]  } 


=  W   (s) 


as  in  equation  (2.14)  above.   For  n=2,  equations  (2.55)  and  (2.56)  agree 
with  equations  (8)  and  (9)  in  Osaki  [34].   Using  the  principle  of 
regeneration,  an  alternative  derivation  of  (2.55)  is  treated  in  the 
Appendix. 


The  3-Out-of-n  System 


For  the  3-out-of-n  svstem,  finding  the  expectation  of  T„ 

J ,  n 

again  necessitates  an  expanded  state  space.   As  before,  let  X(t)  be 
the  state  of  the  system  at  any  moment  t;  e.g.. 
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(i)  X(t)  =  (i,z),  1-  i.  i  i.  ^»  z  >^  0 ,  means  Chat  at  moment 
t,  machine  number  i  is  dovm  and  an  amount  of  time  z,  0  <_  z  <_  t ,  has 
already  been  spent  on  the  repair  of  machine  i. 

(ii)    X(t)  =  (i,j;z),  1  ^  i,  j  1  n,  i  /  j,  0  <_  z  <_  t ,  means 
that  at  moment  t,  machine  number  i  and  machine  number  j  are  dovm,  and 
an  amount  of  time  z,  0  <  z  _^  t ,  has  already  been  spent  on  the  repair 
of  unit  i. 

(iii)   Again,  X(t)  =  (0,0)  means  that  at  moment  t,  no  machines 
are  down. 
For  l-;i,  j<n,  ij^j,  let 


PgCt)  =  P[X(t)  =  (0,0)] 


p^(t,z)cz  =  P[X(t)  =  (i,y)],  z  <  y  <_  z  +  dz 

p^.(t,z)dz  =  P[X(t)  =  (i,j;y)],  z  <  y  £  z  +  dz 

n         / .  V.  n         / .  .  \     T^ 

a=  .1      X.,        X^^^  =   y   A.,    A^^^>  =   V   A  ■ 

j=l  J  j=l  ^             k=l   " 

With  S  =  S(i,j)  =  {i,j:  i,j  =  l,2,...,n;  i  r   j},  the  state  transition 
equations  are 


PQ(t  -f  dt)  =  ?o(t)(l  -  cxdt)  + 

n   r^  ,    .  F.(z  +  dt)  -  F  (z) 

+  I      I   p.(t,z)(l  -  A^^dt)  — — r— -—--— dz  +  o(dt)   (2.57a) 

. :^i  !   1  1  -  F.  Cz; 

1--  ^  1 

(..  1  -  F.(z  +  dt) 

P^(t  +  Qt,z  i-  dt)  =  p^(t,z)(i  -  ),^^dt)  — r^TTz') —  ■""  °<^^^^'  ^1^1^ 

(2.57b) 
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p   (t  +  dc.z  +  dt)  =  p   (t,z)(l  -  A^^J^dt)    ,      \    ,    s + 

iJ  ij  1  -  F.  (z) 

1  -  F.(z  +  dt) 
+  p^(t,z)A^dt  — f-Z-l^— (7) +  o(dt),  i,j  G  S  (2.57c) 

P^(t  +  dt,0)dt  =  pQ(t)A^dt  + 

n   f  ,. .V    F. (z  +  dt)  -  F.(z) 

+  o(dt),  1  1  i  1  n  (2.57d) 

p.^(t  +  dt,0)dt  =  o(dt),  i,j  e  S  (2.57e) 

Pq^O)  =  1  (2.57f) 

Defining  the  functions  {g.(t,z)},  ,  and  q..(t,z),  i,i  z    S,  by  the 

1         1=1         XJ    '    »    »J       I        J 

relationships 

p^(t,z)  =  e      g^(t,z)(l  -  F^(z)),  1  <_  i  _<  n  (2.58) 


p   (t,z)  =  q   (t,z)(l  -  F.(z)),  i,j  €    S  (2.59) 

equations  (2.57a)  -  (2.57f)  may  be  rewritten  using  (2.7)  as 


t 
0       ...."!    -A'"'z 


dp„  n   f    (i) 


dT^  -c»P0<;t)  +  _Z   J   e      g.(t,z)dF.(z),pQ(0)  =  1 


(2.60a) 


0 
+  -~  =0,   1  1  i  .1  n  (2.60b) 


9gi    9gi 


3c    dz 


L't      dz  'ij   '  •     j        -^i 


^^  +  ~-^-  -  -  x'^^J^q.  .(t,z)  +  A.e      ^g  .  ( t ,  z)  ,  i  ,  j  g  S        (2.60c) 


g. (t,Q)  =  p.(t,0) 
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r 
=   Xv(t)    +      I  q..(t,7.)dF.(z),  1  £  i  <  n 

J    n 


j?^i 


(2.60d) 


q.j(t,0)  -  p.j(t,0)  =  0,  i,j  e  S 


(2.60e) 


Now,  for  i,j  (,    S,  put 


Vq(s)  = 


e   PQ(t)dt 


V.(s)  = 


e  ^^(  I  p^(t,z)dz)dt 
0        0 


V,.(s) 


-St, 

e    ( 


p.  .  (t,z)dz)dt 
,   1.1 


G.(s) 


~st  ,   K  , 
e   g^(t)dt 


9.(s)  = 


e  ^*^dF.(t) 


';'.(s)  =  (^.(s  +  >/^^) 


and 


W.   (s)  = 


-•■.r. 


I   V  (s)  +  I      I        V..(s) 


i=0 


i.j  c  S 


Equations  (2.60a)  -  (2.60e)  will  nov  be  used  to  find  W^   (s)  and 

3,n 
~*  -sT 

h,  ^(s)  =  E[c:   '3,n],  as  follows.   Frora  (2.60b) 


g^(t.z)  =  g^(t  -  z),  1  1  i  1  n 


(2.61) 


34 


Using  (2.61)  and  (2.60e)  the  solution  of  (2.60c)  becomes 

q_(t,z)  =  e       g^(t  -  z)  -  e  '    g^(t  -  z)  ,  i,j  e  S     (2.62) 


Substituting  (2.61)  and  (2.62)  into  (2.60d)  and  taking  Laplace 
transforms  of  the  resultant  equations,  obtain 

n 
G.(s)  =  A  V  (s)  +  I      (^  .(s)  -  if-.(s))  •  G.(s),  1  1  i  1  n    (2.63) 

■!=]_        J  J  J 

Transforming  both  sides  of  (2.60a)  with  the  aid  of  (2.61)  yields 


n 


(s  +  cx)Vq(s)  =1+1      ^^(s)    ■    G^(s) 

i=l 


(2.64) 


Equations  (2.63)  and  (2.64)  represent  n+1  relationships  between  the 
n+1  functions  Vq(s),  G^(s),  G^(s),...,G    (s). 

From  (2.58)  and  the  definitions  of  the  V.(s) 


G,  (s) 


V.(s)  =  -i ^  (1  -  .(..(s)),  1  1  i  1  n 


(2.65) 


s  +  A 


and  from  (2.59),  (2.62)  and  the  definitions  of  the  V..(s) 
G.(s)  G.(s) 


s  +  A 


s  +  A 


(2.66) 


Hence,  using  (2.65)  and  (2.66),  for  n  >_  3 


i=l         i,j  E  S   -J 


G.(s)(l  -  ^^..(s)) 


+  a(^^> 


s  +  A 


n  G.(s) 
-  (n-2)j:  -^-^   (l--.(s)) 


i=l  s+A 


(2.67) 
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and 


h*   (s)  =  1  -  sW^   (s)  (2.68) 


E[T^   ]  =  W-   (0)  (2.69) 

3,n     3,n 


where  V  (s)  and  {G.(s)}._^  must  be  found  from  equations  (2.63)  and  (2.64) 


Then 


As  a  check  on  (2.67),  let  a.  =  A  and  F.(-)=F(-),  lli^n. 


)^(s)  =  (^(s),  4^^(s)  =  (^(s  +  (n-l)X),  ^i.(s)  =  4>(s  +  (n--2)A), 


a  =  nX,  A*-^^  =  (n-l)A,  A^^^"*  =  (n-2)A,  i,j  e  S 


For  convenience,  let 

6   -(s)  =  (i)(s  +  (n-l)A) 
n-i 


e   „(s)  =  *(s  +  (n-2)A) 
n-/ 


Y„(s)  =  B^  „(s)  -  3^  ,(s) 
n       n-z       n-i 


In  terms  of  the  above  notation,  equations  (2.63)  and  (2.64)  can  be 
written  in  matrix  notation  as  (equation,  page  36).   Using  induction 
and  Cramer's  rule  on  the  above,  for  n  >_  3 

1  -  (n-l)Y  (s) 

^0^-^^    (s  +  nA)(l  -  (n-l)7  (s))  -  nAB^  ^  (s)  ^^-""^ 

n  n-l 

and  for  1  <_  i  <_  n 

^i  ^ ^ ''  ^  Ts  +   nA)(l  -  (n-l)Y  (s))  "nW'TU)  ^^'^^'^ 

n  n-i. 

iJtev   substituting  (2,70)  and  (2.71)  into  (2.67)  and  simplifying,  it 
■Ls  seen  that 


36 


to 
o 


u 


o 


I 


I 


CO 


I 
C 


I 


ca 
I 


CO 


I 


ca 


+ 
CD 


^ 


?- 


;^ 


CO 


>~ 


CO 


;^ 


?- 


>- 
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as  in  equation  (2.30).   Thus,  for  the  special  case  of  X^  =  X  and 
F.(-)  =  F(-),  equation  (2.69)  agrees  with  equation  (5,19b)  of  Dovmton 

[27]. 

Again,  however,  since  the  state  equations  change  for  each 

value  of  k,  a  general  expression  for 

k,n 

cannot  be  found  with  this  technique.   Nevertheless,  for  any  fixed 
value  of  k,  1  £  k  <^  n,  an  expression  is  obtainable  for  ^^^(s)  and, 

hence,  for 

-sT 
E[e   ''''']  =  1  -  sW^^^(s) 


and 


E[T,   ]  =  W,   (0) 
k,n     k,n 


A  Comparison  Between  the  Semi-Markov  Model  and 
the  Supplementary  Variable  Technique 

General 

Since  both  Downton  [27]  and  Branson  and  Shah  [39]  use  semi- 
Markov  processes  to  model  repairable  systems  with  standbys,  a  comparison 
is  in  order  between  the  semi-Markov  approach  and  the  supplementary 
variable  technique  used  in  this  research. 

The  Approach  of  Branson  and  Shah 

In  the  two-unit  system  considered  by  Branson  and  Shah,  a  unit 
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may  operate  on-line  or  off-line.   The  on-line  unit,  no  matter  which 
unit  it  happens  to  be,  always  fails  with  rate  A   and  is  repaired 
according  to  distribution  function  F  (•)•   The  same  is  true  of  the 
off-line  unit  with  respect  to  X   and  F„  ( • )  .   This  is  not  quite  the 
same  as  the  2-out-of-2  system,  but  it  is  similar. 

By  a  judicious  choice  of  states,  Branson  and  Shah  modeled  their 
two-unit  system  as  a  semi-Markov  process.   Then,  using  a  result 
(Barlow  and  Proschan  [9],  Theorem  2.5,  p.  135)  which  does  not  depend 
on  the  distribution  of  the  time  to  system  failure  itself,  the  MTSF  was 
found,  given  that  the  process  began  in  any  of  the  possible  up-states. 

The  Semi-Markov  Model  of  Downton 

Downton  used  a  semi-Markov  process  to  model  the  k-out-of-n 
system  where  X  and  F(-)  are,  respectively,  the  failure  rate  and  repair 
distribution  for  each  of  the  n  units. 

For  n  >_  k  and  k  =  2,3,  Downton  derived 

-sT^  n 
E[e    *  |0  units  are  down  initially]  (2.72) 

but  these  results  are  special  cases  of  equations  (2.55)  and  (2.68) 

with  X.  =  X  and  F.(-)  =  F(').   In  fairness  to  Downton,  however,  although 

X  X 

(2.72)  is  often  the  case  of  interest,  the  semi -Markov  approach  also 

yields 

-sT, 
E[e   '^^'"ji  units  are  initially  down,  0  <"  i  <  k]  (2.73) 

which  cannot  be  obtained  using  the  supplementary  variable  approach 
above. 
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For  each  value  of  k,  the  appioach  used  by  Dox^7nton  necessitated 

a  matrix  inversion  to  find  the  transform  of  the  appropriate  first 

passage  time  distribution.   Since  the  matrix  to  be  inverted  was 

different  for  each  k,  a  general  expression  for  E[T    ]  was  not  obtained. 

In  the  supplementary  variable  approach,  it  was  the  change  in  the  state 

transition  equations  for  each  k  which  prevented  a  generalization  of 

EFT    ]. 
k,n 

Summarizing  Remarks . 

Because  the  semi-Markov  process  is  defined  on  a  denumerable 
number  of  states,  transforms  of  the  distribution  of  T,  the  time  to 
system  failure,  can  be  found  when  the  process  is  in  any  of  the  up- 
states initially  as  indicated  in  (2.73).   Although  this  is  not  the  case 
with  the  supplementary  variable  approach,  results  implied  by  (2.72) 
can  be  obtained  for  a  more  general  class  of  k-out-of-n  systems  than 
those  treated  by  Downton. 

Since  Branson  and  Shah's  approach  to  finding  the  MTSF  is 
independent  of  the  distribution  of  T  itself,  other  moments  of  the  time 
to  system,  failure  are  not  readily  available.   In  this  sense,  the 
technique  employed  by  Branson  and  Shah  is  less  general  than  that  of 
Downton  or  the  author. 


CHAPTER  3 

SOME  RELIABILITY  CHARACTERISTICS  OF  THE  2-OUT-OF-n 
SYSTEM  AND  THE  TWO-UNIT  STANDBY  REDUNDANT  SYSTEM 


Introduction 

Reliability  attributes  of  the  2-out-of-n  and  two-unit  standby 
redundant  systems  are  treated  in  this  chapter.   Besides  the  distri- 
bution of  the  number  of  repairs  completed  prior  to  system  failure, 
transforms  of  distributions  are  derived  for  the  time  the  repairman  is 
idle  during  T  (the  time  to  system  failure)  and  the  time  spent  on 
repair  during  T.   Assuming  random  failure  and  k-Erlang  repair 
capabilities,  the  distribution  of  T  for  each  system  is  analyzed  and 
the  generating  function  of  the  distribution  of  the  number  of  renewals 
occurring  during  T  is  obtained. 

The  2-Out-of-n  System 
Assumptions  and  Definitions 

The  following  will  be  assumed: 

(1)  The  time  to  failure  of  each  unit  is  independently  and 
exponentially  distributed  with  failure  rate  X  >  0. 

(2)  At  t=0,  all  n  units  are  operating  properly. 

(3)  There  is  only  one  repairman  (with  unlimited  service 
capabilities)  and  the  repair  times  are  independent  random  variables 
which  are  also  independent  of  the  failure  times. 

Define: 

Z(t)     -   the  number  of  units  down  at  moment  t  (Z(t)  =0,  1  or  2) . 

T„      -   the  time  to  system  failure  for  the  2-out-of-n  system.   As 
2  ,n  ^  ^ 
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before,  the  density  of  T„    and  its  Laplace  transform  will 

z ,  n 

be  denoted  by  h„   (t)  and  h„   (s),  respectively. 
2,  n         2,n 

X       -   the  time  elapsed  between  the  completion  of  the  (i-l)-st 

repair  and  the  next  breakdown,  i  ^  1  (the  zeroth  repair  is 

assumed  to  be  completed  at  t=0) .   By  assumptions  (1)  and 

(2),  the  {X.},  i  >_  1,  are  independent  and  identically 

distributed  (i.i.d.)  as  a  random  variable  X  with  density 

function 

-X  X 
fj^(x)  =    ^Q^  ,   X  ^  0 

where  X   =  nX,  n=2,3,.... 


i 


R_,      -   the  time  required  to  perform  the  i-th  repair,  i  ^  1«   By 
assumption  (3),  the  {R.},  i  ^  1,  are  i.i.d.  as  a  random 
variable  with  distribution  function  ?(•)•   Let  (J'Cs)  be 
the  Laplace-Stieltjes  transform  of  F(');  i.e.. 


4.(s)  = 


e  ^'^dF(t) 


0 
Y       -   the  time  to  the  next  failure  when  only  n-1  units  are 

operative.   By  assumptions  (1)  and  (3),  the  density 

function  of  Y  is 

-x^y 

f^Cy)  =  X^e  ■"   ,   y  ^  0 

where  X   =  (n-l)X,  n=2,3 

N       -   the  number  of  repairs  completed  prior  to  system  failure, 

R       -   the  total  time  spent  on  completed  repairs  during  T 
c  z  ,n 


42 


R       -   tha  total  time  spent  on  all  repairs  during  T„ 

i ,  n 

I       -   the  total  time  that  the  repairman  is  idle  during  T„ 

2,n 

A  Conditional  Transform  Approach 

Before  discussing  the  distributions  of  N,  R  and  I,  it  will 
be  instructive  to  consider  a  typical  realization  of  the  process  Z(t) 
(Figure  3.1). 

Note  that 


P[a  repair  interv^al  contains  no  failures] 
=  P[R  <  Y]  = 


e   "■  dF(t)  =  ii\^) 


0 


Clearly,  the  distribution  of  N  is  geometric  and 


P[N  =  i]  =  {C(A^)}^  (1  -  $(X^)),  i=0,l,... 


(3.1) 


Z(t) 


9  .. 


1  -• 


■X. 


"2  ,n 


^3-"^  Y 


Figure  3,1.   A  sample  function  of  Z(t)  in  terms  of  the 
random  variables  X.,  R.,  Y  and  T, 


1'   1' 


2,n- 


Given  that  N  =  i 
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f     R,  +  R^  +  . . .  +  R. ,   i=l,2, . 
12  1 


R  =  < 
c 


^  0,   i=0 


(3.2) 


where  the  {R.},  i  >  1,  are  i.i.d.  as  a  random  variable  R  with 
J     - 


dP[R  <  t]  = 


-At 
e    dF(t) 


and 


(3.3) 


I  =  X,  +  X„  +  ...  +  X. ,, ,   i  >  0 
1    2  1+1     — 


(3.4) 


Hence,  from  (3.1),  (3.2)  and  (3.3) 

N  =  i]P[N  =  i] 


-sR      °°     -sR 
E[e  ''\   =  I  E[e   ^ 


i=0 


i=0  t 


*(s  +  X^) 
*(A^) 


{<^(\^)}\l  -  MA^)) 


1  -  '^(X^) 
1  -  (^(s  +  \^ 


(3.5) 


and  by  (3.4)  and  the  definition  of  the  {X.},  j  2.  •'- 


E[e  ^^\    =  \     E[e  ^^  N  =  i]P[N  =  i] 


i=0 


=  I 


i=0 


s  +  X, 


i+1 


{c^(Ap}"(i  -  -Hx^) 


Xq(1  -  <^(X^)) 
s  +  Aq(1  -  Mx^)) 


(3.6) 


As  a  check,  note  that  (see  Figure  3.1) 


T„    =  Y   +  R  -r  I 
2,n         c 


(3.7) 
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where 


A  e     (1  -  F(y)) 
dP[Y   <  y]  =  — 


1  -  MA^) 


dy,  y  >_  0 


Now,  by  (3.2),  (3.3)  and  (3.4) 


-s^R   -si 
E[e   "  '^  e  ^    ] 


"     -s,R   -s„I 
I      E[e      ^   "  e      ^ 
i=0 


=  I 
i=0 


*(s^  +  A^) 
*(A^) 


N  =  i]'P[N  =  i] 
i+1 


X 


0 


^2  +  ^0 


(*(A^))\l 


(A^)) 


Xq(1  -  (KX^)) 


s^  +  Xq(1  -  4>(s^  +  X^)) 


(3.8) 


Putting  s   =  s„  =  s  in  (3.8) 


-s(R^+I)       Xq(1  -  <^(X^)) 
E[e        ]  = 


s  +  Xq(1  -  ({,(s  +  X^)) 


(3.9) 


By  definition  of  dP[Y   <  t]  above 


-sY*,    [    -st,„.,,*     .     A^(l  -  *(s  +  X^)) 


E[e    l-j   e   dP[Y  i  tl  ■(,+,, (1  -  *(»,)) 
0 


(3.10) 


Corabining  (3.9)  and  (3.10)  with  (3.7) 

-sT  *       -s(R  +1) 

Els     '"^1  =  E[e"®   ]  •  E[e    ^        ] 


XqX^[1  -  Ms  +   X^)] 


(s  +  X^)(s  +  Xq[1  -  ^(s  +  X^)]) 


(3.11) 


and  equation  (3.11)  agrees  with  (2.14). 
Clearly,  the  moments 
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E[{R  }^{I}^].   j,k=0,l,2,..., 
c 

are  available  by  differentiating  (3.8).   Also,  since  R,  the  total  time 
spent  on  all  repairs  before  failure,  satisfies 


R  =  R  +  Y 
c 


by  (3.5)  and  (3.10) 


E[e   ]  = 


+  X. 


(3.12) 


From  (3.12)  it  is  seen  that  R  is  distributed  exponentially  with  mean 
l/X  ,  regardless  of  the  repair  distribution.   For  n=2,  relations  (3.6) 
and  (3.12)  agree  with  equations  (17)  and  (21)  in  Epstein  [31]. 


An  Analysis  of  the  Distribution  of  T^  ^ 

In  this  section,  the  density  of  T^^^.  ^^2,^''^'   ^^'^^  ^^ 
investigated  when  repair  times  follow  the  V-Erlang  distribution  with 

parameter  y  >  0;  i-e.,  when 


*(s)  = 


e  ^^dF(t)  = 


s  +  y 


,  k  >_  1 


(3.13) 


First,  applying  the  shifting  property  of  Laplace  transforms  to  equation 
(3.11)  and  recalling  that  X^  =  (n-i)X,  i-0,1.  obtain 

h    (t)  =  exp{-(n  -  l)Xt}g(t)  (3-1^) 

2,n 


wnere 


g''(s)  =  j   e""'  g(t)dt 
0 

=  h*   (s  -  (n  -  I)X) 
2,n 
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=  n(n  -  1)A^[1  -  (^(s)] 
s(s  +  A  -  nX(};(s)) 

or,  using  (3.13) 


g*(s)  =  C(s)/D(s)  (3.15) 


where 


C(s)  -  n(n  -  l)X^[(s  +  y)^  -  p^]/s 


D(s)  =  (s  +  A)(s  +  y)^  -  nAii'^ 


(3.16) 


and  k>_l,  n^2;  A,y  >  0. 

Note,  however,  that  for  0  <  y  <  A  and  k  >^  l,,n  >  2 


E[repair  time]  =  k/y  >  ((n  -  1)A)  ■""  =  E[Y] 


In  words,  the  expected  repair  time  is  greater  than  the  expected  time 
to  the  next  failure,  given  a  repair  has  just  begun.   Hence,  for 
y  <  A,  the  2-out-of-n  system  will,  "a  priori,"  be  unreliable  and  the 
interesting  case  is  when  y  >_  A  >  0.   Under  this  latter  condition,  the 
following  theorem  will  be  proven: 

Theorem  3.1.   Let  y,  A  >  0  and  let  k  and  n  be  integers  with  k  >_  1 
and  n  >_  2.   If  y  >_  A  >  0,  then  the  k+1  zeros  of  D(s)  =  (s  +  A)  (s  +  y) 
-nAy   are  distinct. 

Proof.   Assume  s  =  r  is  a  root  of  multiplicity  m  >^  2  of  D(s)  =  0; 
then 


^D(s) 

ds-" 


=  0,   j=l,2,...,m;  m>_2 
s=r 


Taking  j=2. 


D'(s)|     =  (r  +  A)k(r  +  y)^  ^  +  (r  +  y)^  =  0 


s-r 
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or 


(r  +  u)^   ^(k(r  +  X)  +  r  +  p)  =  0 


But  r  7^  -p  since  D(-ij)  ^   0,  and  hence 


k(r  +  X)  +  r  +  p  =  0 


or 


r  =  - 


kA  +  y 
k  +  1 


However, 


D(s) 


s  =  - 


kX  +  y 
k  +  1 


=  -{(p  -  X)^^\^  +  nXp^Ck  +  l)^^^}/(k  +  1)^^^ 


<  0,  when  p  >_  A  >  0 

Therefore,  D(s)  cannot  have  any  root  of  multiplicity  m  >_  2  when 

p  >_  A  >  0.  II 

k+1 
Consequently,  if  {r.}._^  are  the  k+1  zeros  of  D(s)  by  using 

(3.14),  (3.15),  Heaviside's  expansion  and  Theorem  3.1 


k+1 


h    (t)  =  exp{-(n-l)At}  •    I      exp(r.t)C(r.)/D'(r.) 


2,n 


j=l 


(3.17) 


where  C(s)  and  D(s)  are  defined  by  (3.16) 


The  Number  of  Renewals  During  T 


2.n 


Suppose  that  at  time  t=0,  a  renewal  process  begins 
which  generates  secondary  failures.   These  secondary  events  are  not 


considered  sei-ious  enough  to  necessitate  inmediate  repair  (Gaver  and 

Luckew  [38])  and  it  is  assumed  that  total  failure  will  not  be  caused 

by  them  alone  (e.g.,  an  oil  leak  or  a  faulty  valve).   However,  when 

failure  of  the  2-out-of-n  system  does  occur,  all  secondary  failures 

are  also  repaired  before  the  system  is  restarted. 

Let  the  independent  times  between  secondary  failures  have  the 

distribution  function  B(")'   Then,  conditional  upon  T    ,  it  follows 

I. ,  n 

from  renewal  theory  that  the  distribution  of  the  number  of  secondary 

failures  during  T    ,  say  M,  is 
z ,  n 


F[M  =  ilT^  n  "  ^^  "  ^i^''-'  "  ^i+l^*"^'  ^  -  °  O.IZ") 


where 


B^^j^(t)  =     B.(t  -  x)dB(x),  i  >_  0  (3.19) 


and  B  (t)  is  the  unit  step  function  at  the  origin. 


Now,  if 


B  (s) 


e  ^'^dB(t) 


0 

and  if  it  can  be  shown  that 


Re(r,)  <  (n  -  1)X,   1  <  j  <_  k+1 


it  follows  that 


exp{-[(n  -  1)>  -  r.]t}B.(t)dt 

[B*((n  -  1)A  -  r  )]^ 

-,j>  — ■  -.  ...1 


(n  -  1)A  -  r. 
J 

for  1  Jl  j  .1  k-M,  i  ^  0.   Using  (3.17),  (3.18)  and  (3.20),  one  may 
calculate 


(3.20) 
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M 


I      z^P[M  =  i] 
i=0 


i=0 


k+1 
=   I 


(B.(t)  -  B.^^(t))h2^^(t)dt 
t:=0 
C(r.)[l-B*((n-l)X  -  r  )] 


j=l  D'(r.)((n-l)A-r.)[l-zB  ((n-l)X-r.) 


(3.21) 


where  0  £  z  j<  1. 

Thus,  in  order  for  (3.21)  to  hold,  there  remains  to  show 

k+1 
Theorem  3.2.   If  p>_A  >  0,  k^l,  n^2  and  if  {r.}.^,   are  the  k+1 

k      k 
zeros  of  D(s)  =  (s  +  A) (s  +  y)   -  nXp  ,  then 


Re(r.)  <  (n  -  1)X,   j=l, 2, . . . ,k+l 


(3.22) 


Proof.   By  Descartes 's Law  of  Signs,  D(s)  has  precisely  one  positive 

(real)  zero,  say  r.   Note  that  D(0)  <  0  and  D((n  -  1)X)  >  0  and 
hence,  it  must  be  that 

0  <  r  <  (n  -  1)X 

and  all  real  zeros  of  D(s)  satisfy  (3.22). 

Now,  with  i  =  /-T,  let  z  =  a  +  3i  be  any  complex  zero  of  D(s) 
such  that  Re(z)  =  a  >^  (n  -  1)X,  and  let 

Q  =  I(z  +  X)(z  +  y)^| 

Since  z  is  a  zero  of  D(s),  it  follows  that 


or, 


(z  +  X) (z  +  p)  I  =  InXy  |  =  nXy 


Q  =  I (a  +  X  +  Bi)(a  +  p  +  6i)  |  =  nXu 
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I     .       .  k|    Ilk 

But  since  | z  z„ |  =  | z^ |  •  | z^ |  and  | z^ |  =  | z^ |   for  any  two  complex 

numbers  z  ,  z  ,  it  follows  that 

I  II  rk 

Q  =  la  +  A  +  6i|  •  |a  +  u  +  Bi| 

=  [(a  +  A)^  +  6^]^^^[(a  +  y)^  +  e^]^/^  (3.23) 


However,  for  a>_  (n  -  1)A,  A,  v  >  0,  k^l,  n>_2  and  for  all  real 
values  of  3, 

2     2  1/2        2     2  1/2 

[(a  +  A)^  +  S  ]  '   >  [(nA)   +  3  ]  '   >  nA  >  0       (3.24) 


and 


[(ct  +  u)^  +  3^]^^^  i  [((n  -  1)A  +  p)^  +  3^]^^^ 

>  ((n  -  1)A  +  y)^  >  u^  >  0       (3.25) 


Combining  (3.24)  and  (3.25)  with  (3.23) 


r/   ,  ,n2  ,  ,2,1/2,/   ,   s2  ,  ^2,k/2    ,  k 
Q  =  [(a  +  A)   +  3  ]    [(a  +  p)   +  3  ]    >  nAu 


whenever  a  >  (n  -  1)A.   Hence,  when  a  =   Re(z)  >_   (n  -  1)A,  the  complex 
number  z  =  a  +  Si  cannot  be  a  zero  of  D(s)  and  the  theorem  is  proven. ] 

From  (3.21),  all  mom.ents  of  M,  the  number  of  secondary 
failures  during  T„   ,  are  available  by  differentiating  G  (z).   In 
particular 

k+1  C(r.)E*((n  -  1)A  -  r.) 

E[M]  =   I  ^ ; J 

j=l   D'(r.)((n-1)A  -  r.)[l  -  B"((n-1)A  -  r.)] 
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The  Two-Unit  Standby  Redundant  System 
Assumptions  and  Definitions 

The  following  will  be  assutaed: 

(1)  The  time  to  failure  of  unit  1  is  a  random  variable  which 
is  independent  of  the  time  to  failure  of  unit  2,  and  switchover  time 

is  instantaneous. 

(2)  At  t=0,  unit  1  is  operating  and  unit  2  is  a  "cold"  back-up. 

(3)  There  is  only  one  repairman  with  unlimited  service 
capabilities.   The  repair  time  for  unit  1  is  a  random  variable  which 
is  independent  of  the  repair  time  for  un: t  2,  and  the  repair  times  and 
failure  times  are  independent  of  each  other. 

Comment:   Once  again,  to  avoid  cum.bersome  or  new  notation, 
the  random  variables  X.,  R. ,  R,.  R,  I,  N  and  M  will  continue  to  be 
used,  even  though  their  meaning  may  be  different  than  in  previous 
sections.   The  reader  is  cautioned  to  make  note  of  this  whenever  these 
quantities  are  defined  below. 

The  following  notation  will  be  used: 

f   -St 
i{'}  -      the  Laplace  transform  operator;  i.e.,  £{f(t)}  =  j   e   f(t)dt. 

0 
Z(t)     -   the  number  of  units  down  at  moment  t  (Z(t)  =0,  1  or  2) . 

T  -  the  time  to  system  failure  for  the  two-unit  standby 
redundant  system.  The  density  of  T  and  its  Laplace 
transform  will  be  denoted  by  h.^,(t)  and  h^^Cs)  ,  respectively. 

i       -   subscript  denoting  unit  i,  i=l,2. 

X       -   time  to  failure  of  unit  i. 
i 
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g.(-) 


probability  density  function  of  X. 


G.(-) 


distribution  function  of  X. 


G.(.) 


1  -  G.(.) 


-   time  needed  to  repair  unit  i. 


f,C) 


F.(.) 


F,(-) 


N 


-  probability  density  function  of  R. . 

-  distribution  function  of  R. . 

1 

=   1  -  F.(-) 

the  number  of  repairs  completed  prior  to  system  failure, 


R 

R. 

1 

I 


the  total  time  spent  on  repair  during  T. 


the  total  time  spent  repairing  unit  i  during  T. 


the  total  time  that  the  repairman  is  idle  during  T. 


A  Conditional  Transform  Approach 

Before  deriving  the  distribution  of  N,  consider  some  typical 
realizations  of  the  process  Z(t)  (Figures  3.2  and  3.3). 
Now,  let 


p  =  P[R^  <  X^]  = 


g2(t)F^(t)dt  = 


f^(t)G2(t)dt 


q  =  P[R2  <  X^]  = 


g^(t)F2(t)dt 


f2(t)G^(t)dt 


and  note  that  (see  Figures  3.2  and  3.3) 


P[N  =  n,  N  odd]  =  p("+l)/2  q(n-l)/2  ^^_^^ ^  n=l,3,5,...  (3.26) 


P[N  =  n,  N  even]  =  (pq)"^^(l  -  p) ,  n=0,2,4,. 


(3.27) 
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z(t) 


2  ■• 


1  ■• 


0 


^1^- x^  —^ x^ ^ x^ ^  h- 


i^— R, 


Figure  3.2 
variable 

are  completed. 


variables  X. , 
1 


A  sample  function  of  Z(t)  in  terms  of  the  random 
R.  and  T,  when  an  odd  number  of  repairs  (3) 


Z(t) 


) 

T^ 

2    ■ 

• 

h 

w   ^         \r           h^ 

1    ■ 
0 

• 

.^     ?  . 

"2 

-^^2    ^ 

^h-* 

^^2- 

— » 

*-\ 

•   •   •  • 

Figure  3.3.   A  sample  function  of  Z(t)  in  terms  of  the  random 


variables  X. 

1 ' 

are  com.pleted. 


R.  and  T,  vhen  an  even  number  of  repairs  (2) 
1 
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or 


-[N  =  n]  =  I  {(-1)^""^  +  l}p^^+l^^'2  q(n--^)/2  (i  _  q)  + 


+  I  {(-1)''  +  l}(pq)"^^l  -  p),   n=0,l,2,...        (3.28) 


As  a  check,  using  (3.28) 


G  (z)  =  I      z"p[N  =  n] 
n=0 


1  -  p  +  p(l  -  q)z 

,     2 
1  -  z  pq 


2        ,   0  £  z  £  1  (3.29) 


and  from  (3.29) 


G^(l)  =  1 


To  find  the  transform  of  the  distrihution  of  T,  a  conditional 
transform  approach  is  used: 

E[e~^  ]  =    y    E[e~^'^|N  =  n,  N  odd]  •  P[N  =  n,  N  odd]  + 
n  odd 

V       -sT  I 
+   2.    ^[^    |N  =  n,  N  even]  •  P[N  =  n,  N  even] 

n  even 

(3.30) 

For  N  =  n  and  N  odd 


(n-l)/2        (n+l)/2 

T  =  X  +      y    X,  .  +      y    x„  .  +  X 


1 


where,  the  {X^  .},  j  >  1,  are  i.i.d.  as  a  random  variable  X,  with 
l.J     -  1 

g,  (t)F,/t) 

dP[X.  <  t]  =  — = dt 

1  -  q 


and   {X      .},    J   1  1.    are   i.i.d.    as   a   random  variable   X^  with 

g2(t)F^(t) 
dP[X„    <    t]    = dt 
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and 


g^(t)F2(t) 


dn\it]  =     ^  _  q 


dt 


Also,    for  N  = 


n   and  N  even 
n/2 


n/2 


T   =   X, 


■"     I        h^-"     ^        ^2,j    ^ 


whare 


*  go(t)F    (t) 

dPLx^  1  t]  =  -  ,.p       dt 


Defining 


gj^(s)    =   Mg^(t)} 


and  combining  (3.26)  and  (3.27)  with  (3.30)  and  the  above  discussion, 
it  follows  that 


h*(s)  =  E[e"^"] 


n  odd 


'^^§2^''^^1^^''^ 


(n+1)/: 


£{g^(t)F2(t)} 


(n-iy2 


■  Ji{g^(t)F2(t)} 


1  -  q 


,(n+l)/2   (n-l)/2(^  _  ^^ 


ft      ^ 
+  g^Cs)    I 


n  even 


'  Mg2(t)F^(t)}  ^ 


n/2 


£{g^(t)F2(t)} 


n/2 


'  £{g2(t)F^(t)}  ~ 


1  -  p 


(pq)''/^l  -  p) 
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^,     £{g2(t)F^(t)}  +  9Ag^(t)F^(t)}    •  ^{g2(t)F^(t)} 

(3.31) 
Equation  (3.31)  agrees  with  equation  (90)  in  Srinivasan  [25]  and 
equation  (14)  in  Osaki  [34]. 

Using  a  similar  approach,  transforms  of  the  distributions  of 
R,  R^,  R  and  I  can  be  found.   In  particular,  for  N  =  n  and  N  odd 

(n+l)/2        (n-l)/2 


j=l    ^'J     j=l    ^'J    ^ 


and  for  N  =  n  and  N  even 


n/2         n/2   ^       ^ 

R  =   I    ^1  i  +  I    ^2  i  +  ^2 
j=l   ^'^        j=l   ^'^    "^ 


where,  the  {R   .},  j  >_  1,  are  i.i.d.  as  a  random  variable  R,  with 

f  (t)G  (t) 

dP[R,  <  t]  = dt 

1  -  p 

and  {R„  .},  j  >_  1,  are  i.i.d.  as  a  random  variable  R  with 
^  >  J  "t 

f,(t)G  (t) 

dP[R„  <  t]  =  ~ i—  dt 

Z  -  q 

and  dP[X^  i.  t] ,  j  =  1,2,  are  as  defined  previously.   Hence,  from 
(3.26),  (3.27)  and  the  above  discussion 


E[e  ^^]  =   I        E[e  ^^|n  =  n,  N  odd]P[N  =  n,  N  odd]  + 
n  odd 


V       — sR  I 
+   2,     E[e    |N  =  n,  N  even]P[N  =  n,  N  even] 

n  even 
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£{g,(t)F_(t)}  •  ii{f.(t)G„(t)}  +  £{g„(t)F  (t)} 
=  L £ -J: 1 _£ ± (3.32) 

1  -  £{f^(t)G2(t)}  •  J2.{f2(t)G^(t)} 


Similarly, 
-sR 


E[e    h    =    I 

n  odd 


£{f^(t)G2(t)} 


(n+l)/2 


p(n+l)/2  q(n-l)/2^^  _  ^^^  ^ 


+   I 


n  even 


£{f^(t)G2(t)} 


n/2 


Jt{g2(t)F^(t)} 
1  -  P 


(pq)"/^l  -  p) 


(1  -  q)£{f^(t)G2(t)}  +  £{g2(t)F^(t)} 
1  -  q£^^f^(t)G2(t)} 


and 


-sR„ 

E[e    ]  =   A 


n  even 


£{f2(t)G^(t)} 


n/2 


(pq)^/^l  -  p)  + 


n  odd 


£{f2(t)G^(t)} 


(n-l)/2 


'  £{g^(t)F2(t)}  ' 


1  -  q 


.   (n+l)/2   (n-l)/2,.     . 
•  p        q        (1  -  q; 


p  £{g^(t)F2(t)}  +  (1  -  p) 
1  -  p  £{f2(t)G^(t)} 


A  check  reveals  that 


-sR, 


E[e  ="]     =  E[e   "] 
s=0 


-sR  -sR„ 

^^     =  E[e   n 


s=0 


=  1 


s=0 


Using  the  above  approach,  it  is  easy  to  calculate 


58 


-s,R.   -s  R  -s  R   -s  R 

E[e      e     ]    =        I        E[e      ^        e  |n  =  n,  N  odd]P[N  =  n,  N  odd]  + 

n  odd 

+   I  E[e      e      |n  =  n,  N  even]P[N  =  n,  N  even] 

n  even 


and  show  that 


-s  R   -s  R 
E[e   -^  ^  e  ^    n 


-s(R+R  ) 
=  E[e    ^   ^  j  =  E[e  ^^] 


s^=S2=s 


as  expected,  since  R  =  R  +  R  . 

The  same  type  of  conditional  transform  approach  is  used  to 

find  the  transform  of  the  distribution  of  I,  the  idle  time  of  the 
repairman  during  T. 

Referring  to  Figures  3.2  and  3.3,  it  is  clear  that  some 

conditional  probability  density  functions  (CPDF's)  must  be  calculated; 

namely,  the  CPDF  of  X  -  R  given  that  X.  >  R. ,  say  g^  „  .    ^  (•), 

^1        A--R^  |X_>R 


and  the  CPDF  of  X^  -  R^  given  that  X^  >  R^,  say  g^  _^    .   ^^ 


(•) 


Using  the  transformation  of  variables  technique  and  the  definitions 
of  p  and  q,  it  can  be  shown  that 


2y  -R  Iy  >r  ^^'^    "^     f T  (x  -  y)g„(x)dx/p,  0  <  y  <  « 
2   1 '  2   1      ^ 


x=y 


and 


%,-R^|x  >R.  '•^■^ 


1   2 '  1   2 


^2^^  "  y)8-j^(x)dx/q,  0  <  y 


x=y 


Letting 


and 


gJl(s) 


y=0 


e  ^^'(J    f^Cx  -  y)g^(x)dx)dy 
x=y 


g*2(s)  =  [   e  ^^(1    f,(x  -  y)g.,(x)dx)dy 


J 
y=0 


and  remembering  that 


it  follows  that 


g^(s)  =  £{g^(y)} 


3S 


Eie    ] 


gl(s) 


y 

n  odd 


§21(3) 


(n+l)/2  f      * 


S^^is) 


(n-l)/2 


•  p(n+l)/2  ^(n-l)/2  ^^  _  ^^  ^ 


11  even 


g2i(s) 


n/2 


g;L2^^^ 


n/2 


(pq)''^^(l  -  p) 


g^(s)[(l  -  q)g2-^(s)  +  1  -  p] 
_      _ 

1  -  g-L2(s)g2j^(s) 


(3.33) 


Since  g^(0)  =  1,  g2i(0)  =  P  and  g3^2^°^  "  '^'  ^'^°"'  (3.33) 


E[e-^] 


=  1 


s=0 


as  it  should  be. 


— s (R+I) 

Although  E[e       ]  is  not  readily  available,  using  (3.31), 

(3.32)  and  (3.33),  by  algebra 


E[T]  =  -  ^   E[e-'^] 


s-0 


d   J.  ,  -sR,  ^  „,  -si,  -, 
=  -  —  {E[e    ]  +  r.[e    } 

QS 


s=0 


•=  E[R]  +  E[T] 
as  expected,  since  T  =  R  +  I. 
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-sT      -sR      -sR^ 
Although  inverse  transforms  of  E[e    ],  E[e    ],  E[e     ], 


-sR. 


2         -si 

E[e     ]  and  E[e    ]  are  difficult  to  calculate  for  particular  failure 

and  repair  distributions,  moments  of  T,  R,  R  ,  R   and  I  are  available 
by  differentiation  and  possibly  with  the  aid  of  numerical  techniques 
(reference  Chapter  4) . 

An  Analvsis  of  the  Distribution  of  T 


when 


and 


In  this  section,  the  density  of  T,  h  (t) ,  will  be  investigated 


-At 


.  (t)  =  Ae  ^  ,   A  >  0,  t  ^  0 


.k-1 


f,  (t)  =  p(yt)    exp(-yt)/(k-l)!,   y  >  0,  k^  1,-  t  >.  0 


for  i  =  1,2.   Defining 


(J)(s)  = 


e  ^'^f.(t)dt,   i=l,2 


it  follows  that 


4>(s)    = 
From    (3.31)    and   the   above 


s  + 


^  ) 


,  k  ^  1 


(3.34) 


h^(s)  = 


A  n  -  6(s  +  A)] 


(s  +  a)   -  A(s  +  a)<^(s  +  a) 


(3.35) 


As  before,  by  the  shifting  property  of  Laplace  transforms  and  equation 
(3.35) 

h^(t)  -  e"^^k(t)  (3.36) 

where 


k  (s) 


e  ^\(t)dt 


0 


^  s   -  Xs(}>(s) 
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or,  using  (3.34) 


k  (s)  =  c(s)/d(s) 


(3.37) 


where 


c(s)  =-  X^(s  +  y)^  -  /} 


k     k 

d(s)  =  s(s  +  \i)      -  Ay 


(3.38) 


The  following  theorem  will  be  useful; 


Theorem  3.3.   If  X,  y  >  0,  the  zeros  of 


k     k 

d(s)  =  3(s  +  y)   -  Xy 


are  distinct  for  k=l,2,.... 

Proof .   The  proof  of  Theorem  3.3  follows  the  same  line  of  reasoning 

as  Theorem  3.1.  1 1 

Thus,  if  {s.l^^j  are  the  k+1  zeros  of  d(s),  then  by  (3.36), 


(3.37),  Heaviside's  expansion  and  Theorem  3.3 


^   k+1  s.t 
h^(t)  =  e"^'  I      e   J  c(s  )/d'(s 
^  -1=1        J      -J 


) 


(3.39) 


v'here  c(s)  and  d(s)  are  defined  by  (3.38). 


The  Number  of  Renewals  During  T 

Again,  properties  of  renewal  theory  may  be  used  to  find  the 
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generating  function  of  the  distribution  of  M,  the  number  of  renewals 
(or  secondary  failures)  during  T. 

As  before,  let  the  independent  times  between  renewals  have  the 
distribution  function  B(').   Then 

P[M  =  iJT  =  t]  =  B^(t)  -  B.^^(t),  i  ^  0  (3.40) 

where  the  B.(')  are  defined  as  in  equation  (3.19)  above. 

It  is  asserted  that: 

k+1  1 

Theorem  3.4.   If  v,  X  >  0,  k_^  1  and  if  {s.}._^  are  the  k+1  zeros   of 

k     k 
d(s)  =  s(s  +  u)   -  Ap  ,  then 

Re(s.)  <  X,   j=l,2, . . . ,k+l 

Proof .   The  proof  is  almost  identical  to  that  of  Theorem  3.2  and  will 
not  be  repeated. 1 1 
Hence,  with 

oo 

B*(s)  =     e"^'^dB(t) 


0 
and  using  (3.39),  (3.40)  and  Theorem  3.4,  one  may  calculate 

oo 
CO  oo        .  . 

G  (z)  =  I      z^P[M  =  i]  =   y   2^  i   (B.(t)  -  B    (t))h  (t)dt 
i=0  i=0    •'     ^       i+i     i 

t=0 

k+1       c(s.)[l  -  B*(A  -  s.)] 

j=l   d'(s.)(>,  -  s.)[l  -  zB*(A  -  s.)] 

•J  J  J 

where  0  <  z  <  1. 


See  Sivazlian  [40]  for  an  analysis  of  a  special  case  of  the 
polynomial  d(s)  which  arises  in  inventory  theory. 


CHAPTER  4 
APPLICATIONS 

Introduction 

In  this  chapter,  applications  of  the  theory  in  Chapters  2  and 
3  are  treated.   In  particular,  numerical  results  are  obtained  for  the 
2-out-of-3  and  3-out-of-4  systems  with  respect  to  an  airport  limousine 
problem  and  some  numerical  techniques  applicable  to  Chapters  2  and  3 
are  mentioned. 

An  Example 

A  city  power  plant  system  or  an  airport  limousine  (also 
referred  to  as  car)  service  are  two  practical  situations  which  can  be 
modeled  as  a  k-out-of-n  system.   For  instance,  if  a  city  has  say  three 
power  plants  of  varying  sizes,  it  is  likely  that  if  any  two  of  them 
are  in  a  failed  state,  the  third  power  station  will  become  overloaded 
and  the  entire  system  will  fail;  i.e.,  a  2-out-of-3  system,  assuming 
one  repair  crew.   A  more  realistic  and  intuitive  example  is  that 
of  an  airport  limousine  service  and,  hence,  discussions  will  be 
confined  to  the  latter. 

Consider  the  manager  of  an  airport  limousine  service  composed 
of  three  limousines.   Assume  that  so  long  as  at  least  two  of  the 
limousines  are  operative,  airport  customers  will  be  inclined  to  use 
the  limousine  service.   However,  if  at  any  time  only  one  limousine  is 


■'"According  to  Professor  0.  I.  Elgerd,  Department  of  Electrical 
Engineering,  University  of  Florida,  the  k-out-of-n  system  is  too 
simplified  a  model  for  power  plant  systems  in  general;  however,  the 
above  discussion  does  provide  the  reader  with  a  feeling  for  what 
k-out-of-n  systems  are. 
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operative,  customers  will  lose  patience  and  choose  other  means  of 
transportation  (e.g.,  busses  or  taxis).   Thus,  the  manager  would  be 
interested  in  knowing  the  mean  time  to  failure  (one  limousine  working) 
for  his  2-out-of-3  system;  i.e.,  E[T2  ^]    in  the  notation  of  Chapter  2. 
In  addition,  the  manager  might  also  be  interested  in  how  much  an 
extra  limousine  would  be  worth  to  him.   In  other  words,  how  does 
E[T    ]  compare  with  E[f   , ] . 

Assuming  limousine  i  has  an  expected  time  to  failure  of 
1/X  weeks  (1  week  =  7  days),  and  the  mean  time  to  repair  limousine  i 
is  p.  weeks,  1  £  i  £  4,  calculations  of 

^'2, 3^°)  =  ^^^2,3^     ^^^     ^,4^°^  =  Ef^3.4^ 

were  made  using  equation  (2.56)  (for  W^  3(0))  and  equations  (2.63), 

(2.64),  (2.67)  and  (2.69)  (for  W.  ,(0)).   The  results  for  several  sets 

3,4 

of  values  of  {X^}  and  {p.}  are  shown  in  Tables  1  and  2  for  three 
situations,  when  all  repair  times  are  (1)  deterministic,  (2)  2-Erlang 
and  (3)  exponential. 

As  indicated  in  Table  1,  for  the  2-out-of-3  system  exponential 
repair  appears  to  be  better  than  deterministic  repair  capabilities 
when  the  mean  repair  times  are  the  sam.e.   This  peculiarity  of  the  2-out- 
of-3  system  was  also  mentioned  by  Downton  [27].   However,  the  standard 
deviation  of  T^  ^  is  slightly  larger  for  exponential  repair  than  it 
is  for  either  deterministic  or  2-Erlang  repair  (see  Table  3). 

A  comparison  of  Tables  1  and  2  v/ould  aid  the  limousine  manager 
in  deciding  the  value  of  an  extra  limousine.   Comparing  cases  1  and  5 
for  exam.ple,  the  manager  presently  (Table  1)  has  three  cars  vjith  mean 


TABLE  1 
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E[T,  J  FOR  DETEEMTNISTIC.  2-ERLANG  AND  EXPONENTIAL  REPAIR 


Case 
No. 


Failure   Rates   and 

Mean  Repair   Times 

for  Unit _i,    i=l,2,3 


Et'^2^ 


"1 


X^  =   1/8,    A^  =   X„   -   1/A 

P3  =   2/7 


y^   -    1/7,    y^ 


. 


X,    =    X 


1/8,    X-   =   1/A 


2        -'-'      3 


p^   =   V2   ^-^   1/^'    ^3   ==   2/'^ 


Deterministic 
B£gair____, 


X„   =    X, 
2  -> 


y^  =  ^2  =  ^3 


1/8 
1/7 


^1  =   ^'2   =   ^3  =   ^/^ 


19.39 


36.45 


80.00 


12.01 


2-Erl3ng 
_j!.epair_ 


19.80 


36.96 


80.67 


12. 3A 


Exponential 
_Re£air__ 


20.  ly 


37.46  ! 


81.33 


12.66 


TABLE   2 

E[f         1    FOR  DETER>IIiaSTIC.    2-ERLA.NG  Als^D   EXPONENTIAL  PJiPAIR 
^  ,  4 — — — 


T 


Case 
No. 


Failure  Rates  and 
Mean  Repair  Times 
,': o -£ JJnit_i_._  i=1.2,3,4 


A_^  =  1/8,  X2=X3-A^=l/4 
u^  =  1/7,  V'2=V3=y4=2/7 


1^4^ 


Deterministic 
Repair 


6   j  X^  =  1/8,  X2=X3=X^=l/4 


y^  =   1/7,  y2=y3=2/7. 


X^=X„-X3=l/8,  X^  =  1/2 

yi=V2=^y3=i/7,  v^  =  1 


h   =  x^  =  X3  =  x^  =  1/^ 


y. 


y-.  =  y/ 


2/7 


118.38 


32. 


49.53 


62.33 


2-Erlang 
Repair 


87,77 


28.41 


40.18 


47.58 


Exponential 
Repair 


72.61 


26.32 


35.68 


40.33 
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TABLE  3 


THE  STANDARD  DEVIATION  OF  T    FOR  THE  FOUR  CASES 
CONSIDERED  IN  TABLE  1 


Standard 

Deviation  of 

T 
2,3 

Deterministic 

2-Erlang 

Exponential 

Case  No. 

Repair 

Repair 

Repair 

1 

19.28 

19.64 

19.98 

2 

36.35 

36.81 

37.27 

3 

79.93 

80.57 

81.20 

4 

11.89 

12.17 

12.45 

times  to  failure  of  8  weeks,  4  weeks  and  4  weeks  and  mean  repair  times 
are  1  day,  2  days  and  2  days,  respectively.   For  discussion  purposes, 
a  "new  car"  has  8  weeks  and  1  day  as  its  mean  failure  and  repair  time, 
and  an  "old  car"  has  4  weeks  and  2  days  as  its  mean  failure  and  repair 
time.   Thus,  if  the  manager  adds  an  old  car  to  his  original  fleet  of 
one  new  and  two  old  cars,  he  gains  with  respect  to  the  MTSF  almost  100 
weeks  under  deterministic  repair,  about  67  weeks  under  2-Erlang  repair 
and  50  weeks  under  exponential  repair.   If,  however,  after  purchasing 
an  extra  old  car,  it  suddenly  starts  taking  a  full  week  (on  the  average) 
to  repair  it  (e.g.,  parts  must  be  ordered),  from  Table  2,  case  6,  the 
MTSF  shows  a  marked  decrease,  but  is  still  better  than  the  MTSF  for 
his  original  fleet. 

Even  if  the  manager  originally  has  a  fleet  of  three  old  cars 
(Table  1,  case  4),  the  purchase  of  an  extra  old  car  (Table  2,  case  8) 
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adds  40  weeks  under  deterministic  repair,  35  weeks  under  2-Erlang 
repair  and  28  weeks  under  exponential  repair,  to  his  initial  MTSF. 

Although  a  more  complete  study  would  be  necessary  to  draw 
legitimate  conclusions,  it  is  clear  that  the  addition  of  one  extra 
limousine  of  relatively  good  quality  substantially  increases  the 
difference  between  E[T   ,]  and  E[T2  ^]. 

Using  (2.56),  the  benefit  gained  by  having  a  repairman  can 

also  be  seen.   In  particular,  if  t    is  the  time  to  system  failure 

z  ,n 

for  the  2-out-of-n  system  without  repair,  as  (}) .  ->  0  it  follows  from 
(2.56) 


r  n 


2,n 


i--.l  ^ 


-1 


n  X. 


rJ.I^J'"^^ 


;vhen  n=3,  values  of  E[t    ]  for  the  four  cases  in  Table  1  are  shown 
in  Table  4.   A  comparison  of  ELT^  ^]    from  Table  1  and  E[t2  ^]    from 
Table  4  illustrates  the  benefit  gained  by  having  a  repairman. 


T/iBLE  4 


E[t„  ^]  FOR  THE  FOUR  CASES  CONSIDERED  IN  TABLE  1 


Case   No. 

E[i^^^] 

1 
2 
3 
4 

4.13 
5.83 
6.66 
3.33 
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Numerical  Methods  Applicable  to  Chapters  2  ar^d  3 

Except  for  particular  cases  (e.g.,  ^^7hen  all  units  fail 
exponentially  and  repair  time  distributions  are  of  the  Erlang  or 
gamma  type),  the  inverse  transforms  of  many  of  the  results  in  Chapters 
2  and  3  cannot  be  found  by  elementary  methods.   The  reader  is  referred 
to  Linton  [41]  for  a  numerical  approach  which  can  be  used  to  plot  the 
density  of  a  random  variable  which  is  characterized  by  an  irrational 
Laplace  transform. 

However,  as  long  as  the  failure  and  repair  distributions  are 
expressable  in  closed  form,  moments  of  the  random  variables  treated 
above  are  available.   For  example,  the  integrals  necessary  to  find 
E[T]  from  (3.31), 


g2(t)F^(t)dt     and     |   g^(t)F2(t)dt 
0  0 

can  in  some  cases  be  computed  directly  using  Gauss-Laguerre  quadrature 

fcrmulas.   In  any  event,  one  may  also  use  a  transformation  like 

X  =  e   to  transform  (0,  «=)  to  a  finite  interval,  and  then  apply 

Simpson's  rule  or  the  Trapezoidal  rule. 

Hence,  although  moments  of  the  random  variables  discussed  in 

Chapters  2  and  3  are  obtainable,  finding  the  distribution  of  a  random 

variable  (or  a  plot  of  its  density)  from  an  irrational  Laplace 

transform  will  depend  on  the  structure  of  the  case  in  question. 


CHAPTER  5 
CONCLUSIONS  AND  AREAS  FOR  FUTURE  RESEARCH 

Coiiclusions 

The  research  reported  in  this  dissertation  has  broadened  both 
the  class  of  systems  and  system  characteristics  which  can  be  treated 
by  mathematical  reliability  theory.   Although  theoretical  aspects  were 
emphasized,  applications  were  also  considered. 

An  investigation  of  methods  for  analyzing  systems  with  standbys 
resulted  in  generalized  theory  for  finding  the  transform  of  the 
distribution  of  T,   ,  the  time  to  system  failure  for  the  k-out-of-n 
system.   The  main  contribution  of  the  research  was  the  computation  of 
the  transform  of  the  distribution  of  T^^^  via  the  supplementary  variable 
technique,  under  a  more  general  set  of  assumptions  than  has  been 
c=nf:id£red  to  date.   Using  the  principle  of  regeneration,  an  alternative 
derivation  of  the  transform  of  the  distribution  of  T^^^  was  also 

obtained. 

To  investigate  system  characteristics  which  have  been  generally 
neglected  in  the  literature,  a  conditional  transform  approach  ■.vas  applied 
to  Che  twc-uait  standby  redundant  system  and  the  2-out-of-n  system. 
Transforms  of  distribution  were  derived  for  the  time  spent  on  repair, 
the  free  time  of  the  repairman  and  the  number  of  repairs  completed.   For 
z\:e    case  of  exponential  failure  and  Erlangian  repair  capabilities,  the 
time  to  system  failure  of  each  system  was  analyzed  and  the  generating 
function  for  the  number  cf  renewals  occurring  during  the  life  of  each 

system  was  found. 

In  the  context  of  an  airport  limousine  problem,  numerical 
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results  V7ere.  calculated  for  the  2-out-of-3  and  3-out-of-4  systems.   The 
value  of  an  extra  limousine  as  well  as  the  benefit  gained  by  having  a 
repairman  were  discussed.   Numerical  methods  applicable  to  finding 
moments  of  random  variables  characterized  by  complex  integrals  were 
also  mentioned. 

Areas  for  Future  Research 
The  k- quj: -of -n  System 

In  Chapters  2  and  3,  whenever  a  k-out-of-n  system  was  treated, 
it  was  always  assumed  that  failures  were  random  and  repair  was  general, 
but  only  one  repairman  was  permitted.   Hence,  two  obvious  extensions 
are  to  allow  both  failures  and  repairs  to  follow  general  distributions 
and  to  allow  several  repairmen. 

For  the  case  of  several  repairmen  and  general  repair  distri- 
butions, the  ability  to  analyze  the  system  will  depend,  as  always,  on 
how  the  states  are  defined.   If  the  supplementary  variable  approach 
were  used,  two  supplementary  variables  would  most  likely  be  required. 

When  general  failure  and  repair  times  are  assumed,  the  problem 
of  finding  the  distribution  of  T  (or  its  Laplace  transform)  for 
even  the  2~out-of-2  system  has  not  been  solved.   One  approach  in 
particular  is  worth  mentioning  because  of  the  insight  it  provides.   For 
the  2-ou!:-of-2  system  with  one  repairman,  let  Z(t)  be  the  number  of 
units  down  (0,  1  or  2)  at  moment  t.   For  reasons  Xi7hich  will  become 
clear  in  a  moment,  assume  that  each  unit  is  repaired  by  a  different 
repairman,  and  for  i  =  1,2,  define 
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fO,  if  unit  i  is  working  at  moment  t 
1.  if  unit  i  is  under  repair  at  moment  t 


Now,  if  T  is  the  time  to  system  failure,  then 

T  =  min{t;Y^(t)  =  Y^Ct)  =  1}  (5.1) 

A  typical  realization  of  the  processes  Y  (t),  Y^Ct)  and  Z(t)  (see 
Figure  5.1)  helps  validate  the  above  relationship. 

Although  meaningful  results  were  not  obtained  in  this  investiga- 
tion, equation  (5.1),  or  some  equivalent  form,  could  perhaps  be  used 
to  find  the  transform  of  the  distribution  of  T  for  this  general 
2-out-of-2  system. 

The  2 -out -of-n  System  and  the  Two-Unit  Standby  Redundant  System 

In  Chapter  3,  the  transform  of  the  total  time  spent  on  repair 
(R)  was  calculated  for  both  the  2-out-of-n  system  and  the  two-unit 
standby  redundant  system  using  a  conditional  transform  approach.   From 
Figures  3.1  and  3.2,  however,  it  is  clear  that  computation  of  the 
appropriate  stochastic  integrals  (Parzen  [42])  would  also  yield  the 
distribution  of  R  for  each  system.   In  particular,  letting 

2.n 
dI(t,y)=dP[!    Z(u)du  _!  ylT    =  t] 

y       y  J  '^j" 

u=0 

and 


d  I(t,y)  =  d  P[    Z(u)du  <_   y|T  =  t] 
u=0 


it  follows  that  for  tha  2~out-of-n  system 
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Y^(t) 


1  •• 


41 

T      - 

<~  y  -^ 

,   ,              Y 

^ 

^1 

^1 

^"l    -" 

—  R^-* 

-R^-> 

^h    -••• 

Y^CO 


1  ■- 


z(t) 

2   + 


R  ■^—  •  •  * 


Figure  5.1.   A  typical  realization  of  the  processes  Y  (t),  Y  (t) 
and  Z(t)  in  terms  of  the  random  variables  X.  (the  time  to 


failure  of  unit  i) ,  R.  (the  repair  time  for  unit  i)  and  T. 


e-^^C 


— sR 

d  I(t,y)h„   (t)dt)dy  =  E[e   ] 


y=0       t=0 
and  for  the  two-unit  standby  redundant  system 
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e-^^( 


sR. 


d  I(t,y)h^(t)dt)dy  =  E[e   ] 


y=0 


t=0 


though  d  Kt.y)  and  d  i(t,y)  are  not  easily  obtained. 


Hence,  even  thoug..  „    -,^,      ^ 
properties  of  the  random  variables 


Z(u)du 


and 


Z (u) du 


0  0 

are  available  by  means  of  Ele   J. 

Although  a  powerful  tool,  the  application  of  stochastic 
integrals  to  reliability  models  has  not  as  yet  been  researched. 

Priority  Models 

A  third  major  research  area  concerns  the  use  of  priorities. 
In  this  dissertation,  it  was  always  assumed  that  units  were  repaired 
in  the  order  in  which  they  failed  (i.e.,  "first-come,  first-served"). 
By  classifying  each  unit  with  a  priority  index,  the  preemptive-resume 
or  head-of-the-line  disciplines  could  be  imposed  and  new  results  might 
be  obtained  (see  Jaiswal  [22]  and  Natarajan  [28],  [29]). 


APPENDIX 


AN  APPLICATION  OF  THE  PRINCIPLE  OF  REGENERATION 

In  the  spirit  of  Muth  [35],  the  principle  of  regeneration 
will  be  used  to  find  the  transform  of  the  distribution  of  the  time 
to  system  failure  for  the  2-out-of-n  system. 

Define 

T„      -  The  time  to  system  failure  for  the  2-out-of-n  system. 
2,n 

-■k 

As  before,  let  h„   (t)  and  h„   (s)  denote  the  density  of 
<i  ,n         '^  >n 

T    and  its  Laplace  transform,  respectively. 

U.      -  The  time  to  failure  of  unit  i,  1  <  i  <  n.   The  density 

-A  .u. 
function  of  U.  is  assumed  to  be  A. e     ,u.  >0,l<i<n. 

Y.      -  The  time  to  the  next  failure  when  unit  i  is  under  repair. 
The  density  function  of  Y.  is  X^"^   e    ^^  »  7^  ^  0> 
1  <^  i  <_  n,  where 

r  •\  ^ 

X^^^    =  I      X. 

X       -  The  time  elapsed  until  the  first  failure  occurs.   The 

""CtX 

density  function  of  X  is  ae    ,  x  >  0,  where 


ij 


a  -      I      X 


3=1  ' 


D       =   T-    -  X 
2,n 


R,  .     -  The  time  needed  to  perform  the  j-th  repair  of  unit  i, 


l£i£n,  j=l,2,....   The  {R.  . } ,  j  >_  1,  are  assumed  to  be 
i.i.d.  as  a  random  variable  with  distribution  function 
F.  (•) ,  1  i.  i  £  n. 
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z(t) 


The  number  of  units  down  at  moment  t  (Z(t)  =0,  1  or  2) . 


Assume  that  the  failure  times  are  independent  random  variables 
which  are  also  independent  of  the  repair  times,  and  all  units  are 
operative  at  t=0. 

Unless  otherwise  stated,  density  functions  will  be  denoted 

* 

by  a  lower  case  f  and  Laplace  transforms  of  densities  by  f  ,  each  with 

an  appropriate  subscript  for  the  random  variable.   The  convolution 

operator  will  be  denoted  by  an  asterisk  (*) . 

~* 

Before  deriving  h„   (s),  consider  a  typical  sample  function 
z  ,n 

of  ZCt)  (Figure  A.l).   Since  T^    =  X  +  D, 

2,n 


^2,n(^)  =  h^'^'^B^'^ 


(A.l) 


and  note  that  for  1  <  i  <  n 


?[U.  <  U^,...,U.  <  U._^.  U.  <  U.^^....,U.  <  UJ  =  ^   (A.2) 


;(t) 


2   .. 


-»-<- 


•2,n 


Figure  A.l.   A  sample  function  of  Z(t)  in  terms  of  the 

random  variables  X,  D  and  t. 

2,n 


Now,  from  the  principle  of  regeneration,  (A. 2)  and  the  laws  of 
probability, 
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.„(.)  =  l^ 


X^^h-^        \l   -   F.(t))  ^ 


+ 


+ 


x=0 


i=l 


h„   (t  -  x) 
2,n 


(A. 3) 


From  (A.l)  and  the  definition  of  X 


(A. 4) 


Letting 


4>^(s)  = 


-St 


(i) 


dF.(t),  ^As)    =  <}).(s  +  X'  O,  1  1  i  1  n 


and  transforming  both  sides  of  (A. 3),  it  follows  that 


(=)  =  i  I  ^ 


n   X.X^^^l  -  li'.(s)) 


'D      a. 


i=l      s  +  X 


(i) 


—  + 


\.^^^   -.   j,  ^^i 


(s) 


(A. 5) 


From  (A. 4)  and  (A. 5) 


2,n 


n   X^X^^-'CI  -  -^.(s)) 

I     ~                  (i)   ~ 
i=l     s  +  X^  ^ 


n 

s  +  y  A.(i  -  v.(s)) 
i=l  ^ 


(A. 6) 


and  equation  (A. 6)  agrees  with  equation  (2.55) 
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